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1.  INTRODUCTION 


System-generated  electromagnetic  pulse  (SGEM  P)  effects  on  space  vehicles  have 
been  addressed  with  extensive  analytical  and  experimental  efforts.  Complex  analytical 
models  have  been  developed  and  tested  against  experimental  results,  and  key  SGEMP 
parameters  and  hardening  methods  have  been  identified.  Considerable  additional  effort  is 
still  needed  in  applying  the  science  to  military  systems,  however,  and  a  means  of  effi¬ 
ciently  identifying  system  weaknesses  to  SGEMP  effects  would  be  valuable.  Present 
analyses  of  systems  are  generally  conducted  by  a  team  of  specialists  employing  analytical 
and  computer-aided  methods  in  conjunction  with  experimental  efforts.  However,  exten¬ 
sive  training  in  the  field  is  required  to  perform  an  analysis,  and  in  fact,  military  systems 
may  evolve  faster  than  the  present  methods  would  permit  estimations  of  hardness  to 
SG  E  M  P. 

This  report  presents  an  SGE  M  P  analysis  system  designed  to  aid  in  rapidly  scoping 
SGEMP  effects.  A  set  of  ana lyt ical  models  is  solved  on  a  programmable  calculator,  giving 
numerical  evaluations  of  important  parameters  that  characterize  SGEMP  response.  The 
user  must  be  familiar  only  with  parameters  describing  the  x-ray  environment  and  system 
geometry.  Details  of  SGEMP  modeling  are  handled  by  the  calculator  codes.  Output 
parameters  such  as  electromagnetic  fields  and  cable  current  driver  values  permit  the  ana¬ 
lyst  to  determine  magnitudes  of  effects  and  thei r  sensitivity  to  various  parameters.  The 
calculator  codes  are  useful  to  the  experienced  SGEMP  specialist,  and  also  enable  the 
novice  to  identify  potential  areas  of  SGEMP  vulnerability  for  a  particular  system  and 
x-ray  environment. 

The  degree  of  accuracy  which  may  be  achieved  on  the  calculator  is  best  illus¬ 
trated  by  comparing  results  with  more  complex  calculations.  Figure  1  shows  peak  external 
replacement  current  magnitudes  on  a  cylinder  subjected  to  a  variety  of  x-ray  spectra  and 
fluence  levels.  The  curves  were  obtained  at  considerable  expense  (-10  hours  of  CDC 
CYBER-176  computer  time  using  300,000  60-bit  words  of  memory),  using  the  2-1/2- 
dimensional  ABORC  computer  code,  and  at  considerably  less  expense  (~20  minutes  of  exe¬ 
cution  t ime  using  120  words  of  memory)  with  the  Texas  Inst ruments  T 1-59  calculator.  The 
agreement  between  the  two  methods  is  a  factor  of  2  over  problem  conditions  ranging  from 
non-space-charge-limited  to  highly  space-charge-limited.  This  quality  of  results  can 
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Figure  1.  Peak  surface  currents  versus  spectrum  temperature  for 
different  fluence  products  in  cal/m.  Solutions  apply  for 
a  pulse  t  ime  At  =  0.54(2irR/c),  where  c  is  the  speed  of  light. 
The  x-ray  pulse  rise  t  ime  was  =  0.6  x  At. 


probably  be  achieved  in  most  SGEMP  areas  discussed  here,  although  certain  calculational 
areas  require  some  additional  development  to  achieve  thequality  shown  in  Figure  1.  Addi¬ 
tional  comparisons  of  the  external  response  of  the  cylinder  with  ABORC  results  are  found 
in  Section  5,  along  with  an  explanation  of  the  scaled  format  in  which  the  results  are 
presented. 

The  analysis  tools  presented  here  form  a  unified  package  which  allows  SGEMP 
assessment  beginning  with  the  x-ray  environment  and  ending  with  cable  current  drivers. 
The  package  is  complete  in  that  each  major  SGEMP  area  can  be  scoped  with  the  codes. 
The  required  analysis  leading  to  estimates  for  fields  and  cable  currents  is  summarized  as 
follows.  Fi rst,  the  x-ray  spectrum  is  attenuated  to  the  satellite  location  of  interest.  Data 
bases  are  used  to  calculate  values  for  electron  photoemission  yield,  average  velocity,  and 
the  dipole  moment  in  dielectrics.  These  results  are  then  employed  in  equations  for  the 
peak  time  rates  of  change  of  the  electric  field  and  the  space-charge  dipole  moment  near 
an  emitting  surface.  Finally,  simplified  expressions  are  evaluated  for  the  cable  current 
driver  terms. 

Emphasis  in  the  present  program  has  been  on  providing  an  analytical  tool  based 
on  the  development  of  the  field  of  SGEMP  analysis  to  date.  The  analytical  models  dis¬ 
cussed  in  this  report  allow  estimates  to  be  made  for  electromagnetic  fields  in  the  non¬ 
linear  fluence  regimes  —  for  example,  where  the  physics  is  understood  but  the  formulas 
are  usually  too  complex  to  be  solved  by  hand  in  a  practical  manner.  Calculator  results  can 
be  obtained  within  a  matter  of  minutes,  with  only  modest  effort  required  of  the  analyst. 
Quick-response  SGEMP  assessments  with  nontrivial  input  parameters  can  be  performed 
with  the  system.  Potential  applications  include  evaluation  of  candidate  x-ray  simulator 
spectra  for  SGEMP  effects  and  the  determination  of  principal  cable  current  driver  mech¬ 
anisms  for  a  given  threat  environment  and  satellite  orientation. 

The  calculator  codes  here  are  a  generalized  extension  of  work  reported  earlier 
(Ref.  1).  The  analysis  system  described  in  Reference  1  was  limited  to  preselected  satellite 
positions  and  emission  materials,  and  stopped  short  of  the  electromagnetic  field  and  cable 
coupling  capablities  available  with  the  present  system. 

Model  development  to  date  has  been  for  solution  of  satellite  materials,  environ¬ 
ments,  and  x-ray  threat  parameters.  The  particular  codes  discussed  here  were  thus 
designed  for  the  parameter  ranges  and  physics  of  the  satellite  problem.  They  can  be 
applied  to  other  problems  as  well,  such  as  missile  and  g round  systems.  However,  some 
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important  physical  effects  will  not  be  calculated  in  those  cases.  For  example,  air  chem¬ 
istry  effects  important  at  lower  altitudes  are  not  included  in  the  models,  and  x-ray  ener¬ 
gies  are  limited  to  1  MeV  in  the  data  base  (which  is  slightly  lower  than  the  range  of 
interest  for  ground  systems). 

The  remainder  of  this  report  is  organized  as  follows.  Section  2  provides  a  review 
of  SGEMP  phenomenology  and  satellite  system  analysis  methods.  Section  3  discusses  the 
approach  to  simplification  of  SGEMP  for  the  calculator  for  each  analysis  step  and  also 
notes  the  related  phenomena  not  covered  by  the  programs.  Section  4  then  gives  an  over¬ 
view  of  the  modeling  approach  for  each  category  of  SGEMP  analysis.  In  Section  5,  capa¬ 
bilities  of  the  resultant  calculator  system  are  demonstrated  through  comparisons  with 
more  complex  calculations.  A  summary  of  capabilities  is  given  in  Section  6,  which 
includes  quantitat  ive  accuracy  est  imates  and  a  list  of  modeling  limitations. 

Details  of  the  modeling  of  effects  are  found  in  Appendix  A.  Complete  user 
instructions  for  the  codes  are  given  in  Appendix  B,  along  with  a  detailed  example  problem. 
Graphs  of  the  photoemission  excitation  parameter  data  base  are  provided  in  Appendix  C. 
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2.  REVIEW  OF  SCEMP  PHENOMENOLOGY  AND 
SYSTEM  ANALYSIS  METHODS 


The  SGEM  P  problem  is  illustrated  in  Figure  2.  An  x-ray  envi  ronment  specified  by 
a  time-dependent  energy  spectrum  and  flux  is  incident  on  a  satellite  system.  Photoemis¬ 
sion  due  to  the  x  rays  occurs  on  and  throughout  the  structure,  inducing  currents  in  the 
electronics.  Satellite  locations  representative  of  important  sources  of  cable  currents  in 
the  system  are  identified  in  the  figure. 

Lack  of  attenuation  on  the  outside  of  the  vehicle  can  cause  very  large  currents, 
which  must  be  prevented  from  coupling  to  the  inside.  Moderate  attenuation  into  large 
cavities  can  cause  substantial  fields  inside  the  vehicle  which  may  drive  currents  on  cable 
shields.  H  ighly  attenuated  x  rays  can  cause  response  i  n  sens  it  ive  locations  such  as  an  elec- 
tronic  equipment  box.  Also,  x  rays  drive  currents  directly  in  the  cables  through  photo¬ 
emission  from  shield  materials,  penetrating  dielectrics,  or  gaps  within  the  shield.  Thus, 
both  high  and  low  x-ray  fluence  locations  can  be  significant  in  causing  currents  in  the  sen- 
sit  ive  electronics. 

Coupling  of  currents  into  cables  occurs  by  processes  of  knock  -on  current,  BA 
voltage,  CV  current,  and  photon-driven  current.  In  this  discussion,  the  word  ’cable"  can 
mean  the  conductor  of  an  unshielded  cable,  the  conducting  sheath  of  a  shielded  cable,  or 
the  outer  conducting  shield  around  a  cable  bundle.  "Direct  drve’  is  also  used  to  describe 
the  photon-driven  current. 

The  knock-on  current  is  the  net  flux  of  photoelectrons  captured  or  emitted  by  the 
cable  sheath.  The  BA  voltage  source  is  due  to  a  t  ime-rate-of-change  of  magnetic  flux 
enclosed  by  a  cable  loop.  The  CV  current  source  results  from  the  t ime-rate-of-change  of 
the  voltage  between  a  cable  and  another  conductor,  where  C  is  the  mutual  capacitance. 
Finally,  the  photon-driven  current  is  the  result  of  the  displacement  of  charge  emitted 
from  cable  conductors  and  lodging  in  dielectric  insulators  (usually  identified  with  the 
inside  of  a  shielded  cable). 

Other  sources  of  currents  caused  by  the  weapon  include  dielectric  charging  and 
penetration  of  charge  due  to  the  enhanced  electron  environment. 

The  contributions  of  the  different  phenomenology  processes  to  the  total  system 
response  are  il lustrated  in  Figure  3.  The  boxes  represent  calculation  steps  which  must  be 
performed  in  analyzing  a  system.  X-ray  transport  and  photoelectron  emission  are  funda¬ 
mental  to  each  cable  driver  process  shown  in  the  second  layer  of  boxes.  Resultant 
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Figure  3.  Phenomenology  processes  cont  ributing  to  system  response 


energetic  electrons  couple  into  cables  and  electronics  by  means  of  electromagnetic  fields, 
direct  x-ray  effects  inside  cable  shields,  and  effects  inside  equipment  boxes  themselves. 
The  total  system  response  is  obtained  by  simultaneously  combining  each  of  the  effects  into 
an  electrical  model  of  the  whole. 
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3.  SIMPLIFICATION  OF  SCEMP  FOR  THE  CALCULATOR 


The  approach  taken  in  model  development  for  the  calculator  is  typical  of  that 
followed  when  a  full-scale  computer  is  being  used.  The  problem  complexity  or  detail  of 
solution  is  adjusted  to  available  resources,  of  which  memory  si/e  and  execution  time  are 
usually  the  most  critical.  Data  bases  are  often  employed  to  describe  essential  inputs 
which  are  not  the  object  of  the  invest  iRation,  and  detailed  solutions  of  models  containing 
the  parameters  under  investigation  are  solved  using  a  combination  of  analytical  and 
numerical  techniques. 

T  he  calculator  has  many  of  the  attributes  of  a  large  computer,  the  principal  dif¬ 
ferences  being  memory  size  and  execution  speed.  In  general,  greater  reliance  is  placed  on 
analytical  assumptions  and  interpolations  than  would  be  necessary  on  a  large  computer, 
but  the  models  can  achieve  considerable  complexity  andstill  be  solvable  within  a  minute 
or  so  of  execution  time.  The  models  described  here  were  developed  employing  key  SC  EM  P 
parameters  for  simplified  conf  igurations,  with  reliance  on  curve-fitting  to  precomputed 
results  limited  to  the  photoelectron  excitation  parameters. 

X-ray  transport  is  computed  for  incident  snectra  using  exponential  attenuation, 
with  a  simplified  expression  for  the  mass  absorption  coefficients  for  all  materials.  Photo¬ 
electron  parameters  are  obtained  from  convolution  of  attenuated  spectra  with  a  precom¬ 
puted  data  base.  Electromagnetic  response  parameters  obtained  from  simplifying  assump¬ 
tions  for  plane  geometry  are  then  applied  to  finite  objects.  Cable  coupling  terms  are 
treated  separately  for  each  effect,  ignoring  synergistic  possibilities.  Knock-on  current  is 
computed  disregarding  the  the  cable's  effect  on  space  charge,  for  example.  CV  and  ita 
coupling  models  employ  a  rod-over-ground  plane  ideal  geometry,  and  the  direct-drive 
model  assumes  a  coaxial  geometry. 

Effects  of  spacecraft  charging  dup  to  the  enhanced  nuclear  elect  ron  en  vi  ronment 
are  not  modeled  for  the  calculator  at  present,  nor  are  discharges. 

The  present  codes  enable  the  user  to  compute  individual  cable  current  drivers  as 
the  highest  level  of  SC,  EMP  analysis.  The  models  stop  short  of  combining  the  drivers  into 
an  electrical  model  of  the  entire  system,  as  illustrated  by  the  lowest  box  in  T  igure  3.  Such 
a  model  would  suffer  greatly  from  oversimplification  with  the  present  calculator  resources 
due  to  the  degree  of  coupling  between  spacecraft  components.  The  cable  current  drivers 
are  considered  to  be  sufficient  here  for  analyzing  SC  E  M  P  effects  with  the  present  limita¬ 
tions  of  calculator  resources. 
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4.  OVERVIEW  OF  SG  E  M  P  MODELING  FO» 
THE  PROGRAMMABLE  CALCULATOR 


This  section  provides  an  overview  of  the  modeling  of  the  different  SG  E  M  P  analy¬ 
sis  steps  for  the  calculator.  A  brief  description  of  calculator  capabilities  is  given  to 
acquaint  the  reader  with  the  size  of  the  programs  used  to  solve  the  models.  The  idealized 
analysis  conf  iguration  is  described,  and  the  modeling  approach  for  each  division  of  SGEMP 
effects  is  highlighted.  Various  comments  on  the  relevance  of  the  analysis  configurations 
and  model  completeness  and  accuracy  are  also  made.  Details  of  the  modeling  are  dis¬ 
cussed  in  Appendix  A. 

4.1  ALLOWABLE  COMPLEXITY  PERMITTED  BY  CALCULATOR  MEMORY 

The  approach  to  the  modeling  has  been  to  restrict  allowable  complexity  so  that 
all  programs  fit  within  the  TI-59  calculator  in  the  standard  off-the-shelf  configuration. 
Additional  memory  obtainable  through  use  of  a  'chip'  is  not  required,  for  example,  and 
neither  is  a  printer.  Extension  of  the  programs  could  be  accomplished  with  more  memory, 
of  course,  and  increased  accuracy  and  user  convenience  would  be  realized.  At  present, 
however,  the  available  memory  is  assumed  to  be  limited  to  the  standard  960  program 
steps.  (A  program  step  represents  a  single  key  stroke  on  the  calculator  such  as  a  multiply 
or  storage  operation.) 

Calculator  memory  can  be  partitioned  variably  between  storage  of  data  or 
program  steps,  with  a  single  data  number  occupying  the  space  of  eight  program  steps. 
Divisions  of  memory  have  been  varied  for  different  aspects  of  the  SGEMP  problem  in  the 
programs  reported  here.  For  example,  photon  attenuation  and  excitation  parameter  calcu¬ 
lations  require  manipulation  of  a  photon  spectrum,  and  the  calculator  stores  considerable 
data  but  employs  relatively  simple  programs.  The  memory  division  is  240  program  steps 
and  90  data  locations  in  those  cases.  Electromagnetic  response  parameters  require  rela¬ 
tively  few  data  inputs,  but  more  complex  programs.  In  that  case,  the  calculator  memory 
is  divided  into  560  program  steps  and  50  data  locations. 

In  all  cases,  storage  is  severely  limited  when  50  to  90  data  locations  are  com¬ 
pared  against  the  currently  available  420,000  60-bit  words  on  the  DNA  CYBtR-176 
computer,  where  mos  t  of  the  world's  SG  EMP  calc  ul  at  ions  have  been  conducted  to  date. 

Considering  the  limitations  of  calculator  memory,  the  SGEMP  problem  has  been 
broken  into  four  divisions  which  are  treated  by  separate  programs.  The  divisions  are: 
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photon  attenuation,  photoelectric  excitation,  electromagnetic  response,  and  cable  coup¬ 
ling.  The  programs  are  executed  sequentially  by  reading  in  a  different  card(s)  between 
divisions.  This  procedure  is  a  tradeoff  between  problem  resolution  and  operational  con¬ 
venience.  Additionally,  a  user  will,  in  many  cases,  be  concerned  with  certain  of  the  prob¬ 
lem  divisions  more  than  others  in  an  analysis.  He  may  attenuate  his  spectrum  through  a 
thermal  blanket  only  once,  for  example,  and  then  evaluate  the  impact  of  changing  the 
emission  material  or  the  incident  fluence  on  fields  in  a  cavity  by  using  the  attenuated 
spectrum  repeatedly. 

4.2  analysis  configuration 

The  idealized  analysis  configuration  is  illustrated  in  Figures  4  and  5.  The  first 
figure  shows  the  parameters  required  to  define  the  system  for  SGEMP  effects.  They 
describe  the  x-ray  environment  and  the  system  geometry  and  material  composition.  The 
second  figure  gives  the  desi  red  outputs  from  the  mode  Is  by  which  SG  E  M  P  effects  are  mea¬ 
sured.  The  four  cable  current  drivers  are  identified  in  Figure  5,  along  with  other  calcu¬ 
lator  outputs  considered  significant  for  determination  of  SG  E  M  P  effects.  Notice  that  out¬ 
puts  are  available  both  exterior  and  interior  to  the  satellite  model. 

The  cylindrical  geometry  shown  in  Figures  4  and  5  is  employed  to  il  lust  rat  e  the 
relevant  parameters  for  SGEMP  effects.  The  calculator  codes  are  not  designed  to  solve 
the  cylinder  per  se,  but  the  results  can,  in  general,  be  applied  with  reasonable  accuracy  to 
that  geometry.  It  is  convenient  that  major  SGEMP  drivers  are  not  strongly  dependent  on 
object  or  cavity  shapes  for  relatively  simple  geometries  (Ref.  2).  One-dimensional,  planar 
geometry  results  often  can  be  modi fied  slightly  and  applied  to  much  more  complex  shapes, 
especialy  near  photoemitting  surfaces. 

4.3  SUMMARIES  OF  MODELS  FOR  INDIVIDUAL  ANALYSIS  STEPS 

4.3.1  X-Ray  Attenuation 

X-ray  attenuation  is  treated  the  same  as  in  most  other  SGEM  P  codes,  except  that 
it  is  scaled  to  the  calculator  memory  size.  A  universal  expression  for  the  absorption 
coefficient  is  evaluated  at  each  photon  energy  in  the  incident  spectrum.  The  coefficients 
of  the  expression  can  be  varied  for  accurately  fitting  data  for  a  given  material.  Values 
employed  here  were  based  on  the  data  of  biggs  and  Lighthill  (Ref.  3).  Accuracy  is,  in  gen¬ 
eral,  within  a  factor  of  2  of  more  detailed  calculations  for  typical  satellite  materials  and 
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CYLINDRICAL  OBJECT, 


Figure  5.  SGEM  P  response  characterization  parameters  from  calculator  programs 
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x-ray  spectra  (see  Section  5).  Exponential  attenuation  through  parallel  plates  is  performed 
based  on  the  thickness  and  density  of  the  material. 

4.3.2  Photoemission  Excitation 

Photoemission  excitation  parameters  (yield,  average  velocity,  and  dipole  moment 
in  dielectrics)  are  obtained  by  convolving  the  attenuated  photon  spectrum  with  energy- 
dependent  data  for  monoenergetic  photons  obtained  from  the  QUICK  E2  photoemission 
code  (Ref.  4).  Data  for  a  given  parameter  for  a  material  resides  on  one  side  of  a  magnetic 
card  and  occupies  one-fourth  the  calculator  memory  when  it  is  read  in.  The  convolution 
program  performs  power-law  fits  to  the  data  base,  giving  excellent  accuracy  over  the 
entire  range  of  1  to  1000  keV.  Data  are  provided  for  representatives  of  the  full  range  of 
satellite  materials  (listed  in  Appendix  C). 

4.3.3  Electromagnetic  Response 

Electromagnetic  response  parameters  are  obtained  for  a  trapezoidal  pulse  of 
emission  current  from  a  single,  isolated,  conducting,  floating,  infinite  plate.  Simplifying 
assumptions  of  monoenergetic  electrons  emitted  with  a  cosine  angular  distribution  which 

are  not  allowed  to  mix  spatially  above  the  plate  permit  expressions  for  maximum  values  of 

•  • 

the  time-rates^of-change  of  the  electric  field  and  space-charge  dipole  moment  (E  and  P)  to 
be  solved  analytically.  (The  approach  is  similar  to  those  discussed  in  Refs.  5  and  6.) 
Identification  of  P  as  the  principal  driver  of  replacement  currents  on  an  idealized  spher¬ 
ical  geometry  leads  to  an  expression  for  the  magnetic  field,  H.  A  quasi-static  assumption 
leads  to  a  peak  value  estimate  of  the  electric  field.  The  resultant  expressions  are  too 
complex  to  permit  practical  evaluation  by  humans,  but  are  well  within  the  calculator's 
capability  to  solve  rapidly.  Comparison  of  results  for  planar  geometry  with  a  more  exact 
particle-following  computer  code  using  actual  electron  spectra  showed  excellent 
agreement. 

4.3.4  Cable  Coupling 

•  • 

Cable  coupling  terms  are  provided  for  knock-on,  BA,  CV,  and  direct-drive 

effects.  Knock-on  current  is  taken  to  be  simply  the  emitted  current  striking  a  cable  shield 

from  a  nearby  wall  minus  the  current  emitted  from  the  shield.  It  is  simplified  by  ignoring 

the  effect  of  the  cable  on  electron  trajectories,  so  the  results  are  most  accurate  at  low 
•  • 
fluence  levels.  RA  drivers  determine  the  current  through  a  low-impedance  loop  using  B 

obtained  from  P,  the  second  t  ime  derivat  ive  of  the  space-charge  dipole  moment,  and  the 
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•  • 

inductance  of  a  cable  above  a  ground  plane.  CV  coupling  employs  the  electromagnetic  E 

term  combined  with  the  cable  capacitance  above  a  ground  plane.  The  direct-drive  term 
provides  charge  flow  based  on  the  dipole  moment  of  electrons  emitted  from  the  outer 
shield  into  the  dielectric  insulator  of  a  coaxial  cable.  Accuracy  of  the  knock-on  and 
direct-drive  currents  is  most  suspect,  but  all  the  terms  give  reasonable  estimates  of 
effects.  The  electromagnetic  coupling  terms  are  essentially  the  same  as  those  commonly 
employed  in  more  detailed  «?ialyses. 

4.3.5  Summary  of  Modeling 

Table  1  summarizes  the  conditions  for  which  the  various  SGEMP  models  were 
developed.  Table  .!  mmarizes  characteristics  of  individual  output  parameters.  Condi¬ 
tions  and  limiting  assumptions  are  listed  for  each  term,  including  the  assumed  geometry  of 
the  analytical  model  and  the  availability  of  a  nonlinear  (space-charge-limited)  solution. 
Whereas  the  system  covers  many  effects  and  conditions  of  interest,  it  must  be  exercized 
with  some  care  as  to  interpretation  of  results  for  a  given  geometry  and  exposure  environ¬ 
ment.  An  overview  of  the  expected  accuracy  of  the  system  is  presented  in  Section  6. 


Table  2.  Output  Parameter  Characteristics  Summary 


Parameter 

E  xterior 
Solution? 

Interior 

Solution? 

Nonlinear 

Solution? 

Geomet  ry 

Emission  current, 

Yes 

Yes 

N/A 

Plate 

Average  velocity,  <v> 

Yes 

Yes 

N/A 

Plate 

Dipole  moment  rate,  P 

Yes 

Yes 

Yes 

Floating  plate 

Electric  field  rate,  E 

Yes 

Large  cavity 

Yes 

Floating  plate 

E  lectric  f  i eld  E 

Yes 

Large  cavity 

Yes 

Floating  plate 

Magnetic  field  rate,  H. 

Yes 

Pillbox 

E  xterior  on  ly 

Sphere  exterior, 
cylinder  interior 

Magnetic  field,  H 

4> 

Yes 

Pillbox 

E  xterior  on  ly 

Sphere  exterior, 
cylinder  interior 

Knock-on  current,  lknoclc 

Yes 

Yes 

No 

C^ble  over 
ground  plane 

CV  current,  1  • 

Yes 

Yes 

E  xterior  on  ly 

Cable  over 
ground  plane 

Current  from  HA  voltage,  I  * 

H 

Yes 

Yes 

E  xterior  on  ly 

Cable  over 
ground  plane 

Direct-drive  current,  lnD 

Yes 

Yes 

N/A 

Coaxial  cable 
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5.  SAM  PLE  SCEM  P  ANAL  YSES:  COMPARISON  OF  CALCULATOR 
AND  MULTIDIM  ENSIONAL  COMPUTER  CODE  RESULTS 


E  xa triple  SG  E  M  P  analysis  steps  performed  with  the  calculator  codes  are  described 
in  this  section,  and  several  comparisons  are  made  with  more  detailed  calculations.  Accu¬ 
racy  of  x-ray  attenuation  through  representative  satellite  materials  is  compared  with 
results  from  QU  1C  KE2.  Exterior  respons e  is  calculat ed  for  a  cylinder  ove r  wide  ranges  of 
x-ray  parameters,  and  results  are  compared  with  the  two-dimensional  ABORC  code  (Ref. 
7).  A  1-m  length=di  ameter  inner  cavity  is  analyzed  and  compared  with  ABORC,  also.  Fin¬ 
ally,  an  electronics  equipment  box  conf  iguration  is  analyzed  on  the  calculator,  giving  rela¬ 
tive  cable  pickup  magnitudes  for  a  typical  x-ray  environment.  The  results  serve  both  as 
illustrations  of  capabilities  and  as  sample  problems  which  aspiring  analysts  can  use  in 
practicing  the  routine. 

5.1  BLACKBOUY  ATTENUATION  BY  SATELLITE  MATERIALS 

Calculations  of  blackbody  spectra  attenuated  through  representative  satellite 
materials  are  compared  with  QUICKE2  code  results  in  Figure  6.  Agreement  between  the 
methods  is  within  a  factor  of  2,  even  for  the  most  difficult  case  of  large  attenuation 
through  high-Z  materials  (shielded  box).  The  calculator  spectra  were  defined  with  14 
energy  bins,  and  the  QU  1C  K  E2  spect  r  a  had  50.  The  hi  n  select  ion  is  much  more  critical  for 
the  highly  attenuated  cases  than  for  the  moderately  attenuated  cases.  Notice  that  a 
multiple-plate  case  was  performed  on  the  calculator,  also. 

5.2  EXTERIOR  RESPONSE  OF  A  CYLINDER 

The  exterior  response  of  a  cylindrical  object  has  been  calculated  with  the  T  1-59 
codes  for  a  wide  variety  of  incident  x-ray  environments.  The  results  have  been  compared 
with  published  solutions  from  the  ABORC  code  over  ranges  of  x-ray  fluence,  t  ime  history, 
and  energy  spectrum  (Ref.  8).  The  calculator  codes  employ  specific  geometry  assumptions 
involving  planar  and  spherical  surfaces,  whereas  the  ABORC  results  are  for  a  cylinder. 
Previous  analytical  studies  showed  only  moderate  geometry  sensitivity  of  SC  F  M  P  effects 
for  relatively  simple  geometries  (Ref.  2).  Comparison  between  the  calculator  'pseudo- 
sphere*  and  the  ABORC  cylinder  should  not  be  affected  strongly  by  geomet  ry  differences. 
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ATTENUATED  FLUENCE/INC I  DENT  FLUENCE 


Figure  6.  X-ray  shinethrough  fraction  calculated  with  TI-59 
attenuation  program,  compa redto  QUICKE2  code  results  for 
representative  satellite  materials 


Comparisons  between  the  calculator  anrl  the  multidimensional  code  results  are 
shown  in  Figures  7  and  8.  The  graphs  are  presented  in  scaled  format  (Kef.  8)involvingthe 
products,  OK  ("fluence  product'),  E  ,K  ('electric  f  ield  product'),  and  normalized  incident 
pulse  time  history,  AT /2n  K/c,  where  AT  is  the  x-ray  pulse  width  and  2n  K/c  is  the  time 
requi  red  for  an  elect  romagnet  ic  wave  tot  ravers  e  thecircumferenceof  the  cylinder. 

Figure  7  shows  the  effects  of  varying  the  incident  pulse  width  for  a  constant 
spectrum  at  different  fluence  levels,  and  Figure  8  shows  the  effects  of  varying  the  inci¬ 
dent  fluence  for  different  x-ray  spectra.  (Figure  1  shows  the  effect  of  varying  the  inci¬ 
dent  spectrum  temperature  at  constant  pulse  width  for  different  fluences.)  This  scaled 
format  extends  the  applicability  of  the  results  to  a  wide  range  of  object  sizes.  As  an 
example  of  the  use  of  the  scaled  quantities  indicated  in  the  figures,  Figure  7  can  be  used 
to  estimate  the  peak  response  of  a  2-m-radius  object  when  excited  by  a  200-nsec  photon 
pulse  at  a  fluence  of  5  x  10  cal/m  .  For  K  =  2  m,  the  5x10  cal/m  isofluence  line 
corresponds  to  the  bottom  curve  in  the  figures.  For  100  nsec  and  K  =  2  m,  the  abscissa  is 
(10  sec)  x  (3  x  10  m/sec)/(2n  x  2  m),  or  2.39.  The  corresponding  total  current  on  the 
object  is  about  0.5  amps.  The  magnetic  field  is  1/2  it  K ,  or  4  x  10  amp/m. 

Agreement  between  the  calculator  and  the  2-1)  code  results  is  seen  to  be  excel¬ 
lent  in  the  space-charge-limited  regimes,  as  well  as  in  the  non-space-charge-limited 
regimes  which  involve  relat  ively  short  pulse  lengths  compa  red  to  object  dimensions.  The 
agreement  begins  to  deteriorate  particularly  for  low-fluence,  long-pulse-length  problems 
(not  ice  the  I  owe  r  right -hand  port  ion  of  Figure  7)  due  to  limit  at  ions  of  the  present  mode  I  i  ng 
when  the  electron  cloud  extent  becomes  comparable  to  object  dimensions.  This  limitation 
is  discussed  in  Keference  9  in  relat  ion  to  large-sca  I  e  calc  ul  at  ions.  By  virtue  of  the  scaled 
format,  the  inaccurate  region  of  the  parameter  range  can  be  seen  to  become  greater  as 
the  object  size  is  decreased. 

Overall  accuracy  of  the  system  can  be  described  as  excellent  for  the  external 
SG  E  M  P  problem,  however.  Errors  remaining  can  be  reduced  s  ignif  icant  ly  with  moderate 
improvements  in  the  modeling. 

5.3  CAVITY  KESPONSE 

The  internal  response  of  a  cylindrical  cavity  has  been  quantified  with  the  calcu¬ 
lator  codes  and  compared  against  ABOKC  results.  The  cavity  is  a  1-m  lengthaliameter 
empty  cylinder  with  3-mil  aluminum  walls,  and  is  representative  of  the  larger  empty 
volumes  found  in  satellites.  Moderate  attenuation  of  x  rays  caused  by  the  thin  walls 
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QUICKE2  EXCITATION 
PARAMETERS 


IELD  =  1.2  x  10'6  C/cal 
vz  =  2.1  x  10'7  m/sec 


RE-03191  2ttR/C 


Figure  7.  Peak  surface  current  versus  characteristic  pulse  time 
for  various  fluence  products  using  a  2-keV  blackbody  spectrum. 
Fluence  products  $R  in  cal/m  are  noted  bv  each  curve. 

The  x-ray  pulse  rise  time  was  -  0.6  x  At. 


results  in  significant  signals  inside.  Moderate  space-charge-limiting  occurs  at  the  assumed 

_  1  n 

fluence  level  of  10  cal/cm  .  The  incident  spect  rum  is  a  5-keV  blackbody  with  rise  t  ime 
and  pulse  width  (FWHM )  of  10  and  17  nsec,  respectively. 

Results  from  the  two  methods  of  calculation  are  compared  in  Figure  9.  The 
n-59  excitation  parameters  (le  and  <v>)  agree  very  well  with  results  from  the  QUICKER 
code,  mainly  because  differences  should  be  due  to  energy  bi  n  I  i  mi  tat  ion  on  ly.  1  he  elect  ro- 
magnetic  field  peak  values  are  less  accurate  due  to  geometry  problems.  The  calculator 
results  are  three  to  five  times  higher  than  the  more  exact  ABORC  values  due  to  present 
limitations  of  the  simple  models  for  this  case.  In  the  test  problem,  the  electrons  can  tra¬ 
verse  the  entire  cavity  during  the  pulse  width  of  the  x  rays,  whereas  the  analytical  models 
are  exact  only  for  cases  in  which  the  electron  cloud  is  small  compared  to  object  dimen¬ 
sions.  Other  geometry  effects  contribute  to  errors,  including  cylinder  aspect  ratio  and 
electron  angular  distributions.  The  agreement  between  the  calculations  is  sufficient  for 
problem  scoping,  however  (notice  the  calculator  gives  upper-bound  values),  and  variations 
of  input  parameters  would  tend  to  give  reasonably  accurate  assessment  of  tradeoffs  in 
designing  for  hardening. 

The  comparison  problem  shown  here  points  out  the  limitations  of  the  present 
solutions  for  small  internal  cavities  and  low  fluence  levels.  Increasing  the  fluence  would 
result  in  improved  correlation  of  the  parameter  values  between  ABORC  and  the  calculator 
by  reducing  the  space-charge  cloud  dimensions  relative  to  the  cavity  size.  Similarly,  a 
tlecreased  x-ray  pulse  width  would  result  in  greater  accuracy. 

5.4  BOX  CURRENTS 

The  calculator  codes  have  been  applied  to  analysis  of  SC  EMI'  effects  on  a  printed 
circuit  (PC)  board  inside  a  typical  satellite  electronics  equipment  box  (see  Figure  10). 
Worst-case  estimates  were  obtained  of  SGEMP  currents  flowing  on  individual  grounded 
lands  on  the  fiberglass  board  for  a  given  x-ray  envi ronment.  Normally,  such  an  analysis 
might  be  conducted  with  simple  analytical  estimates  once  the  photon  attenuation  and 
excitation  parameters  were  established,  with  some  degree  of  accuracy,  by  a  computer  cal¬ 
culation.  The  results  shown  in  the  figure  can  be  obtained  in  about  10  minutes  with  the 
calculator,  with  the  analyst  having  to  press  only  a  few  buttons  rather  than  having  to  know 
many  different  formulas.  Parameter  variations  for  hardening  measures  or  arbitrary  spec- 
t  ra  can  be  conducted  on  the  box  conf  iguration  wit  h  ease,  a  Iso. 
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PARAMETER 

QUICKE2/ 

ABORC 

CALCULATOR 

ELECTRON  VELOCITY, 

v  (m/sec) 

3.3  x  107 

3.3  x  10/ 

ELECTRON  YIELD,  Y  (C/cal) 

1.8  x  10‘7 

2.3  x  10'7 

EMISSION  CURRENT,  Jg 

2 

(amp/m  ) 

100 

135 

ELECTRIC  FIELD  RATE, 

E  (V/m/sec) 

1  2 

3.5  x  10 

1.5  x  1013 

PEAK  ELECTRIC  FIELD, 

E  (V/m) 

3.8  x  104 

1.1  x  1C5 

MAGNETIC  FIELD  RATE, 

H  (amp/m/sec) 

9.2  X  108 

3.8  x  109 

PEAK  MAGNETIC  FIELD, 

H  (amp/m) 

12 

34 

Figure  9.  Comparison  of  calculator  code  results  with 
<711 1C  K  E  2/A  BOR  C  results  for  interna  I  cavi  t  v  ana  lys  is 


and  the  worst-casp  M  A  is  caused  by  illumination  from  the  bottom,  as  shown. 


6.  SUMMARY  OF  CAPABILITIES  ON  THE  CALCULATOR 


A  summary  of  SGEMP  phenomenology  capabilities  on  the  TI-59  is  shown  in  Table 
3.  The  parameter  ranges  for  which  the  models  were  specifically  developed  are  ind'rated. 
The  codes  can  be  applied  outside  the  indicated  ranges  (except  for  the  energy  range),  but 
the  user  must  be  careful  to  interpret  the  results  properly.  For  example,  mass  absorpton  of 
x  rays  due  to  more  than  15  mils  of  tantalum  can  become  very  delicate  with  only  14  energy 
bins,  and  answers  may  be  very  much  dependent  on  the  particular  energy  bins  selected. 
Attempts  to  quantify  overal I  accuracy  of  the  system  for  SGEMP  effects  within  the  desig¬ 
nated  ranges  are  made,  as  shown  in  the  column  labeled  ‘Accuracy.’ 


Table  3.  Summary  of  SG  E  M  P  Phenomenology  Capabilities  onthe  TI-59 
Parameter  Ranges:  Satellite  locat ions 


X- 

X- 

X- 

ray  energies  from  1  to  1000  keV 
ray  f luence  up  to  1 0 _ *  ca  I/cm2 
ray  pulse  widths  from  2  to  200  nsec 

Calculation 

Accuracy 

Principal  Caveats 

X-ray  attenuation 

Within  factor  of 

2  of  OUICKE  for 
SGEMP  spectra 

Energy  bins  must  be  selected  intelligently 

Compton  effects  ignored 

Exercized  only  for  representative  materials 
to  date 

Photoemission 
excitat  ion 
parameters 

Withi n  factor  of 

2  of  QUfCKE 

Energy  bin  selection 

E&M  response 

Within  factor  of 

Low-fluence,  long-pulse-length  geometry 

pa  rameters 

4  of  ABORC 

Calculator  results 
are  upper  bounds 

effects  not  treated  accurately  (includes 
most  box  effects) 

Moderate  space-charge-limiting  and  angular 
effects  in  cavities  not  treated 

Cable  coupling 

Accuracy  not  tested 
but  calculator  tech¬ 
niques  are  commonly 
used  in  analyses  with 
present  degree  of 
accuracy 

Oirect-drive  est  imate  crude 

Knock-on  current  not  self-consistent 
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Notice  from  the  table  that  factor-of-2  accuracy  can  be  achieved  with  intelligent 
energy  bin  selection  for  the  first  two  categories  of  calculations.  This  is  sufficiently 
accurate  for  scoping  of  effects.  Determination  of  elect  romagnetic  response  parameters  is 
a  little  less  certain,  but  the  calculator  gives  upper  bounds  on  results  in  every  case.  Major 
problems  are  in  the  regimes  of  low  fluence  and  small  dimensions  and  can  be  corrected  with 
simple  analytical  estimates  by  the  experienced  analyst. 

Cable  coupling  estimates  are  accurate  for  cable  lengths  short  compared  to  x-ray 
pulse  lengths.  The  direct-drive  model  is  presently  limited  to  a  simple  est imate  based  on 
the  dipole  moment  from  a  planar  surface.  Accuracy  of  this  model  has  not  been  compared 
with  more  detailed  calculations,  but  better  than  an  order  of  magnitude  is  expected  from 
it.  The  knock-on  current  onto  a  cable  shield  does  not  take  into  account  the  effect  of  the 
cable  itself  on  fields  in  a  cavity,  so  this  dr iver  can  be  low  by  large  amounts,  especial  ly  for 
low-impedance  cables  at  high  fluence. 

In  general,  the  accuracy  for  the  f  irst  three  categories  in  Table  3  is  such  that  ana¬ 
lytical  estimates  of  effects  can  be  upgraded  considerably  through  use  of  the  system  com¬ 
pared  to  present  quick-response  analysis  methods.  The  system  handles  the  most  sensitive 
and  time-consuming  calculational  areas  with  the  greatest  accuracy,  and  areas  of  weakness 
can  be  compensated  for  through  relatively  simple  considerations.  For  example,  inaccuracy 
at  low  fluence  inside  a  cavity  will  likely  be  as  great  due  to  complex  cavity  geometry 
effects  not  treated  in  that  case  as  due  to  modeling  assumptions  used  here. 
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DETAILS  OF  SG  t  M  P  MODELING 
FOR  THE  PROGRAMMABLE  CALCULATOR 


This  appendix  contains  details  of  the  modeling  used  in  the  four  divisions  of 
SGEMP  calculations  treated  here.  Section  A.1  describes  the  photon  attenuation,  A. 2  the 
excitation  parameters,  A. 3  the  electromagnetic  response  parameters,  and  A.4  cable  coup¬ 
ling.  The  physics  and  modeling  for  each  category  are  described  and  compared  against 
more  exact  calculations. 

A.1  X-RAY  ATTENUATION 

Analytical  expressions  for  cross  sections  have  previously  been  developed  (Ref. 
A-1)  and  are  the  basis  for  x-ray  attenuat  ions  i  n  standard  photon  energy  deposit  ion/elect  ron 
emission  codes  (Ref.  A-2).  The  expressions  employ  summations  of  inverse  powers  of  pho¬ 
ton  energy  in  the  photoelectric  regime,  where  the  coefficients  are  based  on  least-square 
fits  to  available  experimental  and  detailed  catculat ional  data.  The  present  work  simply 
reduces  the  method  even  further  to  a  general  approach  suitable  to  rapid  and  convenient 
use  on  a  programmable  calculator.  Limited  storage  capability  prevents  expeditious  hand¬ 
ling  of  large  data  bases,  and  some  degree  of  inaccuracy  is  necessarily  introduced.  The 
simplified  method  and  degree  of  accuracy  are  discussed  below. 

An  expression  which  can  reasonably  approximate  the  photoelectric  absorption 
coefficient  in  the  energy  range  f  =  1  to  1000  keV  is 

k  (Y ) 

<(<f  ')  =  ~  (cm2/g)  ,  (A-1) 

Z 


where  Z  is  the  atomic  number  of  the  material  and  <  is  the  universal  function: 

u 


< 

u 


£  RiH(Y~Yoi) 

i=1  £3 


£  =  rf/Z2  ,  y  =  <f yz2-24  . 


(A-2) 


(A-3) 


Bj  is  a  constant  coefficient  for  an  energy  level  (spe  Table  A-1),  H(x)  is  the  step  function, 
and  Yq|  *  rfyj/Z 2,24  where  is  approximately  the  energy  of  the  ith  absorption  edge  of 

a  material. 


Table  A-1.  Universal  Attenuation  Function  Coefficients 
and  Cutoff  Energies 


Edge 

Z  <29 

l  >  29 

Bi 

(cm2/g) 

Y0i  Bi 
(keV«)m2/g) 

Y0i 

(keV) 

k 

i 

8.6 

4.8  x  10‘3 

12.6 

4.5  x  10-3 

J i 

2 

0.82 

4.8  x  10"4 

1.25 

6.7  x  10-4 

m 

3 

0.13 

1.2  x  10-4 

0.14 

1.2  x  10"4 

n 

4 

0.035 

3.3  x  10"5 

0.010 

1.9x  10‘5 

It  is  important  to  remember  that  only  photoelectric  absorption  is  included  in  the 
expression.  AuRer  absorption  (low  energies)  and  Compton  scattering  (high  energies)  are 
not  treated.  These  effects  tend  to  be  unimportant  for  SC  E  M  P  in  the  energy  range  consid¬ 
ered  here  and  so  are  not  included. 

A  comparison  of  the  universal  attenuation  coefficient  with  more  exact  curves  is 
shown  in  Figure  A-1  for  several  elements  with  a  wide  range  of  atomic  numbers.  Greatest 
errors  occur  for  high-Z  materials  (see  gold).  In  general,  the  absorpt ion  is  approximated  to 
within  a  factor  of  2  of  the  more  exact  curves  until  Compton  scattering  becomes  dominant 
(at  high  energies).  The  approximate  treatment  has  not  been  compared  with  the  more 
accurate  curves  for  every  element,  but  accuracies  similar  to  those  exhibited  are  expected 
for  other  cases. 

The  approximate  method  appears  to  give  photon  attenuations  that  are  suitable 
for  the  present  work.  Errors  introduced  will  tend  to  be  smoothed  out  for  most  spect  ra  due 
to  integration  over  energy  for  most  SGEMP  effects.  Caution  is  recommended  in  using 
monoenergetic  spectra  with  this  method,  since  worst-case  errors  would  occur  with  them. 

The  capacity  of  the  simplified  approach  to  handle  all  elements  with  a  single  alge¬ 
braic  formula  and  a  small  data  base  makes  it  ideal  for  use  in  a  programmable  calculator. 
Additionally,  the  coefficients  of  Table  A-1  can  be  changed  easily  for  more  accurate  fits 
for  a  given  element.  Those  shown  are  chosen  to  represent  all  elements  without  being 
changed  between  calculations. 

It  should  be  noted  that  the  agreements  shown  in  Figure  A-1  were  obtained  with 
coefficients  in  the  attenuation  expression  which  are  material-dependent.  At  present,  the 
calculator  program  employs  a  data  card  for  low-Z  (Z  <  29)  materials  and  a  second  card  for 
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Figure  A-1.  X-ray  photoelectric  absorption  cross  sections  obtained 
using  a  single  expression  for  all  elements.  The  energy-absorption 
cross  sections  of  Biggs  and  lighthill  (Ref.  1)  are  shown  for  comparison. 


high-Z  materials.  Original  plans  for  a  single  set  of  coefficients  for  all  materials  were 
implemented,  but  considerably  better  acuracy  was  realized  when  attenuation  data  for 
ranges  of  material  were  fit  more  exactly. 

The  attenuation  program  is  also  capable  of  treating  multiple  plates  through 
repeated  applicat  ion,  employing  the  output  spectrum  from  one  pi  ate  as  the  i  nput  spect  rum 
to  the  next.  Compounds  are  calculated  by  breaking  the  elements  into  separate  plates 
according  to  their  mass  ratios. 

A. 2  EXCITATION  PARAMETERS 

Excitation  parameter  calculations  are  required  for  the  electric  yield,  average 
velocity,  and  dipole  moment  for  electron  spectra  emitted  into  cable  dielectrics.  A  data 
base  for  these  quantities  has  been  developed  for  representatives  of  the  full  range  of  satel¬ 
lite  materials  (carbon  through  gold),  including  a  few  composites  (SiOj,  Mylar,  fiberglass). 
Values  for  the  excitation  parameters  were  calculated  and  tabulated  as  a  function  of  mono- 
energetic  x-ray  energy  using  the  QUICKE2  photoelectron  emission  computer  code.  The 
results  for  x-ray  spectra  are  obtained  by  convolving  the  monoenerget  ic  values  with  prop¬ 
erly  attenuated  photon  spectra  using  appropriate  algorithms  (Ref.  A-3). 

The  accuracy  of  the  data  base  employed  in  the  calculator  codes  is  essentially  the 
same  as  that  of  the  QUIC  «  E2  code  itself.  The  latter  has  been  shown  to  provide  factor- 
of-2  agreement  with  experimental  data  and  detailed  Monte  Carlo  calculations  for  wide 
ranges  of  photon  and  material  parameters  (Ref.  4-4),  and  is  generally  considered  the 
industry  standard  for  SGEMP  analyses.  Some  inaccuracies  in  interpolating  the  data  base 
values  over  the  range  from  1  to  1000  keV  occur  because  the  calculator  programs  permit  a 
maximum  of  14  energy  bins.  The  accuracy  of  the  excitation  parameters  is  usually  within  a 
factor  of  2  of  QUICKE2  calculations  for  the  same  spectrum,  however,  because  power-law 
fits  to  the  data  base  are  employed. 

F  igure  A-2  shows  a  comparison  of  results  for  excitation  parameters  obtained  with 
the  calculator  codes  and  with  QUICKE2  for  blackbody  spect  ra.  The  agreement  is  consid¬ 
erably  better  than  factor-of-2  over  a  wide  energv  range.  Errors  due  to  inexact  curve¬ 
fitting  at  different  energies  may  be  averaged  out  for  a  spectrum  of  x  rays.  I  he  T 1-59 
results  employed  14  energy  bins  for  photons,  while  QUICKE2  results  were  obtained  using 
50  bins.  Appendix  C  contains  details  of  the  data  base  parameters. 
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ELECTROMAGNETIC  RESPONSE  PARAMETERS 


A. 3 

A  simple  analytical  formulation  is  derived  for  estimating  the  maximum  rates  of 

•  • 

change  of  the  space-charge  dipole  moment  (P)  and  the  surface  elect ric  f  ield  (E)  for  an 
SG  E  M  P  excitation  of  an  inf inite  plane  i  n  f  ree  space.  The  quantities  P  and  E  for  a  planar 
geometry  are  useful  parameters  for  comparing  the  relative  SG  E  M  P  exterior  responses  of 
real  systems  due  to  excitation  pulses  with  different  pulse  shapes,  fluences,  and  electron 
emission  yields  and  energies.  For  an  emission  pulse  whose  intensity  increases  linearly 

from  zero  to  a  maximum  at  time  t^  and  is  then  constant  for  the  remaining  time  of  inter- 

•  • 

est,  the  peak  values  of  P  and  E  and  the  times  of  their  occurrence  can  be  calculated  by 
simple  mathematical  expressions  which  can  be  readily  incorporated  into  a  programmable 
calculator. 

An  expression  for  the  peak  electric  field  at  the  surface  of  a  photon-illuminated 

plate  is  also  described.  Although  this  field  is  not  particularly  useful  for  est  imating  SG  E  M  P 

•  • 

responses,  which  are  primarily  functions  of  B  and  E,  it  is  a  convenient  quantity  for  illus¬ 
trating  the  intensity  of  the  photon  pulse. 

A. 3.1  Foreword 

There  are  two  dominant  mechanisms  for  coupling  the  SCE  M  P  fields  into  an  elec¬ 
tronics  system.  One  is  CV,  where  C  is  the  capacitance  between  a  cable  and  a  nearby 

conducting  surface  of  the  system  and  V  is  the  rate  of  change  of  the  voltage  between  the 

•  • 

cable  and  the  walk  For  a  fixed  geometry,  CV  is  proportional  to  E,  where  E  is  the 
average  electric  field  between  the  cable  and  the  wall.  Since  the  largest  E  fields  on  an 
SG  E  M  P-excited  body  usually  occur  just  above  the  illuminated  surface,  a  one-dimensional 
planar  calculation  of  E  adjacent  to  the  illuminated  surface  should  be  a  measure  of  the 
worst  CV  coupling  that  could  result  from  the  photon  pulse  and  target  material  under 
investigation. 

The  other  coupling  mechanism  is  the  B  coupling  into  the  circuit  due  to  the  mag¬ 
netic  fields  associated  with  the  SG  E  M  P  structural  return  currents.  In  Reference  A-5,  it 
was  shown  that  the  structural  return  currents  (and,  hence,  the  surface  H  fields)  at  a  point 
on  the  outside  of  a  body  that  is  not  too  close  to  the  illuminated  area  of  the  body  are 
reasonably  proportional  to  P,  the  rate  of  change  of  the  space-charge  dipole  moment  at 
the  illuminated  surface.  Thus,  the  peak  value  of  P  and  its  time  derivative  P  provide 
useful  measures  of  the  relative  B  coupling  for  a  system  that  is  exposed  to  different 
photon  pulses.  A  minor,  but  sometimes  useful,  additional  application  of  the  time  to  the 
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peak  P  is  a  guide  for  choosing  suitable  time  steps  for  more  detailed  computer 
calculations. 


In  Section  A. 3. 2,  the  basic  model  used  in  this  report  to  estimate  the  peak  P  and 
E  is  described,  and  the  general  mathematical  equations  for  an  arbitrary  emission  current 
time  history  Uem(t)!  are  derived.  These  equations  are  then  integrated  for  a  trapezoidal 
emission  time  history.  When  the  peak  P  occurs  in  the  linearly  increasing  or  constant 
regions  of  the  trapezoidal  pulse,  the  peak  magnitude  of  P  and  the  time  of  its  occurrence 
can  be  obtained  in  simple  mathematical  closed  forms.  Also,  a  simple  algorithm  gives  an 
approximate  value  for  the  peak  E. 

In  Section  A. 3. 3,  a  comparison  is  given  of  results  from  these  simple  formulas  with 
results  from  a  one-dimensional  particle-pushing  code  which  uses  a  distribution  of  electron 
emission  energies. 

A. 3. 2  Description  of  the  Model 

In  a  one-dimensional  planar  problem  with  the  other  boundary  at  infinity  (an  open- 
circuit  diode),  the  electric  field  at  the  emitting  surface  (E)  is  just  1/e  times  the  total 
charge  densty  per  unit  area  (Q)  that  has  been  emitted  from  the  surface  and  has  not  yet 
returned  to  it.  Early  in  the  pulse,  before  any  charge  has  returned  to  the  surface. 


E(t)  = 


em 


(O  dt  . 


( A—  4) 


Thus,  particles  that  are  emitted  at  time  t  initially  feel  this  electric  field.  Moreover,  until 
these  particles  overtake  some  particles  that  were  emitted  at  anearliertime,  or  unt  il  some 
later-emitted  particles  overtake  those  emitted  at  time  t,  it  is  rigorously  correct  that 
these  particles  wil  I  continue  tofeel  thesame  constant  field  E  (t).  Forthet  imes  of  interest 
for  this  analysis,  it  is  assumed  that  all  emitted  particles  continue  to  experience  the  magni¬ 
tude  of  the  E  field  into  which  they  were  first  emitted  [ E (t)],  even  though  some  spatial 
mixing  of  the  particles  does  occur.  (This  is  the  'gravitational'  approximation  used  in  Ref¬ 
erence  A-6  to  develop  a  thin  boundary  layer  SGEMP  driver  for  a  finite-difference  code.) 

The  next  assumption  of  this  model  is  that  all  of  the  emitted  electrons  have  the 
same  initial  velocity  normal  to  the  emission  surface  (v,„  ).  It  turns  out  that  these  two 

ave 

approximations  give  quite  good  results  up  to  the  time  of  the  peak  P  and  fairly  good 
results  up  to  the  peak  E  as  ong  as  the  energy  spectrum  of  the  emitted  electrons  is  not 
too  broad  with  a  long  range  of  emission  velocities.  With  the  above  approximations,  the 
equation  for  P(t)  is 
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P(t)  =  f  [v  -  E(t')  -  (t-t')l  |  (t')  dt'  , 

L  ave  m  J  e m 


(A-5) 


where  e  is  the  absolute  value  of  the  electronic  charge  and  m  is  the  mass  of  the  electron. 
In  Eq.  A-5,  the  term  in  the  large  brackets  is  the  velocity  at  the  time  t  of  the  electrons 
that  were  emitted  at  t  ime  t'  into  the  constant  electric  field  E  ( t  "  ) . 

A. 3. 3  E  xpression for  P 


For  an  emission  current  )  em(t)  that  increases  linearly  from  zero  to  a  valueof  J 


at  t  =  t. 


P'^'m  vave  “ 


1 2  e  l_  t5  \ 

2t  ~  2  I  ‘ 

1  c  m  40  t  / 
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If  P  reaches  a  maximum  before  t^,  the  maximum  will  occur  at 
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(A-7) 


and  the  peak  value  of  P  can  be  obtained  by  substituting  tm1  in  place  of  t  in  Eq.  A-6. 

For  small  values  of  1^,  the  value  of  tm1,  as  calculated  by  Eq.  A-7,  can  be 

•  • 

greater  than  t^,  so  the  maximum  in  P  occurs  after  t^.  The  equation  for  P(t)  during  the 
constant  part  of  the  trapezoidal  emission  pulse  is 


The  terms  in  the  first  square  bracket  are  the  contribution  to  P  from  the  electrons  emit¬ 
ted  before  t^,  and  those  in  the  second  square  bracket  are  the  contribution  from  the  elec- 
t  rons  emitted  after  t^  (during  the  constant  portion  of  the  emission  pulse). 

The  maximum  value  of  Eq.  A-8  occurs  at 


and  the  maximum  value  of  P  is  obtained  by  substituting  tm2  for  t  in  Eq.  A-fl. 


(4-9) 
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For  very  small  values  of  |  with  a  given  pulse  width,  the  peak  value  of  P  will 

occur  after  the  constant  portion  of  the  trapezoidal  pulse.  It  is  possible  to  obtain  an 

« 

expression  for  P  during  the  linearly  decreasing  leg  of  the  trapezoidal  pulse.  However,  the 
resulting  equation  for  the  time  when  P  is  maximum  is  fourth-degree  and  would  have  to  be 
solved  by  iteration.  Even  though  an  iteration  is  relatively  simple  on  a  large  computer,  the 
storage  requirements  would  probably  exceed  the  capabilities  of  present-day  simple  pro- 
grammable  calculators.  Fortunately,  if  the  peak  P  occurs  after  t^,  one  would  make  only 
a  minor  error  if  one  assumed  that  the  emission  pulse  were  constant  out  to  the  halfway 
point  on  the  linearly  decreasing  portion  of  )  and  used  Eqs.  A-9  and  A-8  to  obtain  the 
time  of  occurrence  and  the  magnitude  of  the  peak  P. 


A. 3.4  E  xpression  for  P 

The  s  implest  way  of  getting  P  ( t )  is  just  to  differentiate  P(t).  Fort<t., 


P(t)  =  )  v 
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Equation  A-10  has  a  maximum  at 
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If  this  t  is  less  than  t-j, 


max 


Vaye  /  Vave  'ijY'3 
*1  \  «m  * 


(A-12) 


For  t1  <  t  <  t2. 


P(t)  =  I 


m 


■  )_e  t 


)_e 


ave 


__  m  r  t 


t  t 


em 


i) 


(A-1  3) 


1  he  first  two  terms  in  Eq.  A-13  (the  constant  terms)  are  the  value  of  P{  t  )  at  t  =  t 

last  term  is  always  subtractive  so  the  peak  of  P(t)  occurs  at  t  =  tj. 

Therefore,  for  P  ,  calculate  t  from  Eq.  A-11.  If  this  t  is  less  than  t,, 
'max  ..  1 

tute  into  Eq.  A-10  or  Eq.  A-12  to  get  P  .  If  that  t  is  greater  than  t,,  evaluate 

max  1 

from  Eq.  A-10  with  t  =  t^. 


The 

substi- 

P 

max 
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A. 3.5 


Expression  for  E 

By  differentiating  Eq.  A-4  with  respect  to  time,  one  obtains 


E(t,  .1  lm(.)  ;  (A-14) 

that  is,  E(t)  is  directly  proportional  to  the  emission  current  up  to  the  time  that  any 
emitted  electrons  return  to  the  emission  surface.  If  this  time  for  the  first  electrons  to 
return  is  greater  than  t1  (that  is,  the  time  when  J  first  reaches  its  peak  value  of  J  m), 
the  peak  value  of  E  is  obviously 


E 

m 


c 


m 


(A-15) 


If  the  first  electrons  return  to  the  emitting  surface  before  t^,  the  curve  of  E(t)  will 
start  to  drop  below  the  linear  curve  of  ((1/e)  !em(t)l  and  then  reach  a  maximum  value. 
Unfortunately,  the  present  model  cannot  be  used  toobtain  a  rigorous  value  for  the  peak  E 
because  the  peak  E  occurs  after  considerable  spatial  mixing  has  occurred  between  elec¬ 
trons  emitted  at  different  times.  Consequently,  several  different  algorithms  were  tried 
with  the  present  model  to  find  an  approximate  method  of  estimating  a  conservatively 
large  value  of  the  peak  E. 

A  relatively  simple  approach  that  gives  somewhat  conservative  (large)  values  for 
•  •  • 
the  peak  E  is  to  assume  that  the  peak  E  occurs  when  P  returns  to  zero.  From  Eq.  A-6, 

P  returns  to  zero  at 


Thus,  by  the  present  assumption,  the  peak  E  would  be  given  by 

L,v  =  7  •  «»<*„*>  * 

ma  x  €  €ftf  m3 


( A-1  f») 


(A-17) 


The  main  rationale  for  this  approach  is  that  it  gives  reasonably  good  comparison 
with  results  from  a  one-dimensional  particle-pushing  code,  as  shown  in  Section  A. 3. 7. 
However,  a  plausibility  argument  is  that  the  dipole  moment  P  is  maximum  when  P  =  0. 
Thus,  on  the  average,  the  emitted  electrons  are  at  the  maximum  distance  from  the  emit¬ 
ting  surface.  Thereafter,  the  emitted  particles  will  return  very  rapidly  to  the  emitting 
surface  and  prevent  E  from  increasing  significantly  more. 


The  general  formula  for  the  peak  steady-state  E,  from  Reference  A-7,  is  (in  M  KS 

units): 

E  m  =  /Emr0  ,  (A-18) 

where  F  is  a  dimensionless  number  that  is  a  function  of  the  emission  spectrum  and  angular 
distribution,  m  =  electron  mass,  rQ  =  electron  emission  rate  =  ) e/e,  is  the  characteristic 
velocity  of  the  electron  spectrum,  and  e  is  the  permittivity  of  free  space.  Reference 
A-7  evaluates  F  for  four  rather  different  spectra,  as  summarized  in  Table  A-3. 

Fortunately,  the  F  terms  differ  by,  at  most,  a  factor  of  1.7  for  these  four  energy 
distributions.  Thus,  the  maximum  steady-state  E  fields  should  differ  by  at  most  /l  .7  - 
1.3.  If  one  just  considers  the  cos  6  distributions  (cases  2,  3,  and  4),  the  F  terms  differ  by 
only  1.5  and  the  E  terms  by  /l . 5  =  1.22.  For  average  results,  one  can  use  F  =  3.0.  Also, 
for  a  cos  0  distribution,  v^  =  /2  vnorma|* 


Table  A-3.  Parameters  Describing  Analytical  Electron  Energy  Spectra 
and  Resultant  Space-Charge-L  imiting  of  Electric  Fields 

Case 

1 

2 

3 

4 

Description 

Monoenergetic 
(energy  =  w.,) 
normal  to 
emission  plate 

Monoenergetic 
(energy  =  w1) 
cos  0  angula  r 
distribution 

Linear  t  imes 
exponential 
energy  distri- 
tion,  cos  0 
angular  dis¬ 
tribution 

E  xponential 
energy  dis¬ 
tribution, 
cos  0  angular 
dist  ri  but  ion 

dn/dw  Delta  function  r^5(w-w^ )  exp(-w/w^ )  —  exp(-w/w.j) 

at  w1(=r05(w-w1)|  w!j  W1 


r 


Figure  A-3  shows  calculated  time  histories  of  E  using  a  one-dimensional  particle¬ 
pushing  code  for  an  electron  spectrum  with  an  average  vnorma|  —  3.54  x  107  m/>ec  and  a 
trapezoidal  emission  pulse  (10-nsec  rise  time),  with  four  values  of  peak  emission  rate.  The 
asymptotic  values  of  the  dynamic  calculations  agree  remarkably  well  (perhaps  partly  for¬ 
tuitously)  with  the  analytical  steady-state  values,  considering  that  the  value  of  F  =  3.0 
which  was  used  is  just  an  average.  Anyway,  the  asymptotic  values  scale  as  /j"^,  as  the 
analytical  expression  predicts. 

Several  other  calculations  using  different  emission  spectra  and  emission  rates 
showed  comparable  agreement,  always  using  F  =  3.0. 

In  conclusion,  it  appears  that  one  is  well  justified  in  using  the  steady-state  for¬ 
mula  (Eq.  A-18)  for  the  peak  E  with  an  average  F  =  3.0.  If  the  pulse  length  is  less  than  the 
space-charge  plasma  period,  this  procedure  overestimates  the  peak  E.  For  example,  if  the 
10  amp/m  pulse  had  ended  at  4  nsec,  where  the  dynamic  calculations  stopped,  the  ana¬ 
lytic  formula  would  have  predicted  E  =  /o.4  x  3.1  x  107  =  1.96  x  107  V/m,  compared  to 
the  dynamic  value  of  1.85  x  107  V/m.  Also,  the  expression  does  not  apply  in  the  non- 
space-charge-limited  regime. 

A. 3. 7  Comparison  to  Results  from  1-0  Particle-Pushing  Code 

The  one-dimensional  particle-pushing  code  used  for  comparison  to  the  above 
simple  model  emits  an  arbitrary  (input)  number  of  particles  per  time  step  with  a  distribu¬ 
tion  of  normal  velocities.  The  total  charge  emitted  per  time  step  is  proportional  to  the 
instantaneous  emission  current,  and  is  distributed  among  the  emitted  particles  consistent 
with  an  arbitrary  (input)  electron  emission  energy  spectrum. 

Results  from  this  code  for  various  trapezoidal  emission  pulses  and  peak  emission 
currents  are  compared  in  Table  A-4  to  results  from  the  present  simplified  model.  The 
electron  energy  spectra  used  in  the  particle-pushing  code  covered  the  range  from  zero  up 
to  130  keV,  with  the  maxima  occurring  at  from  1  to  5  keV.  As  can  be  seen  from  Table 
A-4,  the  values  for  the  peak  P  as  calculated  by  the  simplified  model  are  in  quite  good 
agreement  with  the  results  from  the  particle-pushing  code  (within  about  10%). 

It  is  interesting  that  the  first  ten  problems  shown  in  Table  A-4  used  an  early 
version  of  the  particle-pushing  code  which  did  not  distribute  the  emission  particles  below  1 
keV  because  that  was  as  low  as  the  QUIC  K  E2-generated  emission  spect  ra  extended.  Thus, 
the  velocities  of  the  emission  particles  were  grouped  relatively  closely  around  the  average 
electron  velocity.  The  outcome  is  that  the  P  results  from  the  simplified  model  agree 
almost  exactly  with  the  particle-pushing  results,  as  one  would  expect.  For  the  later 
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Figure  A-3.  Electric  field  at  emission  surface;  1-0  dynamic 
particle-pushinp,  calculation  and  steady-state  analytical  results 
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problems  in  Table  A-4,  the  emission  particles  were  distributed  linearly  down  to  zero  veloc¬ 
ity.  Thus,  the  spread  in  particle  velocities  about  the  average  is  greater.  The  result  is  that 
the  simplified  model  gives  consistently  low  values  of  the  peak  P  relative  to  the  values 
from  the  particle-pushing  code,  but  they  are  still  within  about  10%  of  each  other.  More¬ 
over,  for  comparing  relative  effects  of  different  spectra  and  emission  pulses,  the  simpli¬ 
fied  results  would  give  almost  the  same  results  as  the  particle-pushing  code. 

The  peak  E  values,  as  calculated  by  the  simplified  model,  are  consistently  larger 
than  the  values  from  the  particle-pushing  code,  as  intended,  except  for  the  one  case  in 
which  the  peak  E  from  the  particle-pushing  code  occurred  after  t,.  In  that  case,  the 
agreement  is  exact,  as  it  should  be.  The  largest  discrepancy  is  about  41  %,  which  is  not  too 
bad  considering  that  the  peak  E's  cover  a  range  of  over  three  orders  of  magnitude. 
Moreover,  for  relative  comparison  between  spectra,  the  agreement  is  better  because  the 
errors  in  relation  to  the  particle-pushing  code  are  a  Iways  inthesamedi  rection. 

A  subset  of  the  results  of  Table  A-4  is  plotted  in  F  igure  A-4  for  i  I  lust  rat  ion  pur¬ 
poses.  (The  highest  current  data  points  are  for  a  slightly  different  energy  spectrum  than 
the  others,  and  so  the  curve  is  dashed.)  Excellent  agreement  between  the  two  methods  of 
calculation  is  shown  over  more  than  four  orders  of  magnitude  of  emission  current. 

In  summary,  the  simplified  models  provide  a  convenient  method  for  estimating 

•  • 

the  peak  values  of  P  and  E  and  thei  r  t  imes  of  occurrence  for  a  wide  range  of  spect  ra  and 

•  • 

pulse  shapes.  The  absolute  values  of  the  peak  P  and  E  from  the  simplified  model  are 
within  about  10  and  40%,  respectively,  of  the  corresponding  values  from  the  one¬ 
dimensional  particle-pushing  code.  Peak  E  values  are  consistently  on  the  conservative 
(high)  side,  but  peak  P  values  are  generally  slightly  low.  For  relative  comparisons  among 
different  spectra  and  pulse  shapes,  the  agreement  of  the  simple  model  with  the  particle¬ 
pushing  code  is  even  better  than  the  absolute  agreement. 

The  above  equations  for  the  peak  P  and  E  and  their  times  of  occurrence  are 
quit,e  straightforward  for  programming  on  the  T 1-59  calculator,  with  the  only  required 
input  quantities  being  the  maximum  emission  current  density  (J  ),  the  time  to  this  peak 
emission  current  (t ^ ),  and  the  average  velocity  of  the  emitted  electrons  normal  to  the 
emitting  plate.  For  convenience,  this  average  normal  velocity  can  be  related  to  the  aver¬ 
age  energy  of  the  emitted  electrons,  assuming  a  particular  angular  distribution  such  as  a 
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curves  are  based  on  the  10-nsec  rise  t  ime  r 


A. 3.8 


Estimate  of  Absolute  Magnitude  of  Surface  Skin  Currents 


Skin  currents  on  the  exterior  surface  of  an  illuminated  conducting  body  scale 
approximately  linearly  with  the  rate  of  change  of  the  space-charge  dipole  moment,  P,  in 
the  illuminated  region  (Ref.  A-5).  However,  for  application  toquick  assessments  of  radia¬ 
tion  hardness,  it  would  be  more  useful  if  the  value  of  P  could  be  used  to  estimate  the 
absolute  magnitude  of  the  skin  currents  which  it  would  produce  on  a  body. 

Other  researchers  (Refs.  A-8,  A-9)  have  addressed  this  problem  by  considering 
the  skin  currents  induced  on  a  conducting  sphere  of  radius  R  by  a  single  point  charge  q 
moving  radially  outward  from  the  sphere  with  velocity  v.  When  the  point  charge  is  close 
to  the  surface  of  the  sphere,  the  skin  current,  I  0  (amp),  on  the  sphere  that  crosses  a  cone 
with  the  half-angle  0  about  the  veloci ty  vector  v  ( F igure  A-5)  i s 


I, 


“(■ 


cos  0  + 


/2 


vq 


0  \  /l  -  cos  0 /  2R 

At  the  0  =  90°  point,  Eq.  2  reduces  to 


(A -19) 


lQ  =  0.70  7^.  (A-20) 

In  the  more  usual  situation  with  an  illuminated  sphere,  the  emitted  charge  q 
would  be  spread  over  a  finite  area,  probably  a  hemisphere,  rather  than  being  a  single  point 
charge.  At  a  given  point  on  the  sphere  —  e.g.,  the  0  =  90°  point,  one  would  expect  the 
skin  current  to  be  somewhat  larger  when  the  same  amount  of  emitted  charge  is  spread 
over  a  finite  area  because  part  of  the  emission  a  rea  would  be  closer  to  the  point  where  the 
skin  current  is  monitored.  To  estimate  the  magnitude  of  this  effect,  a  calculation  has 
been  made  for  the  skin  current  crossing  a  cone  with  an  angle  0  when  the  charge  q  is 
emitted  radially  outward  with  velocity  v  with  a  uniform  emission  intensity  o  over  a  spher¬ 
ical  end  cap  that  subtends  a  half-angle  a  at  the  center  of  the  sphere  (Figure  A-6).  The 

m 

density  o  is  given  by 


a  =  — — — q- - .  (A-21) 

2itR  (1  -cos  u  ) 
m 

The  skin  current  is  obtained  by  calculating  the  time  rate  of  change  of  the  surface  charge 
on  the  sphere  inside  the  cone  angle  0. 


I 


I 

Li 
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If  a  and  <}>  are  the  polar  and  azimuthal  coordinates  of  a  point  on  the  sphere  in 
the  emitting  area  (Figure  A-6)  and  the  charge  q  is  at  a  radius  b  slightly  greater  than  R,  the 
equation  for  the  surface  current  is 


( A-2  2) 


C  =  s  in  as  in  0'  cos  <j>  ♦  cos  a  cos  6 


( A — 2  3) 


and 

d2  =  R2  ♦  b2  -  2b  RG  .  (A-24) 

Equation  A-22  has  been  integrated  numerically  for  two  different  values  of  b 
slightly  greater  than  R  since  the  integrand  becomes  infinite  for  some  angles  when  b 
exactly  equals  R.  The  integrated  results  for  the  two  values  of  b  were  close  enough 
together  to  give  confidence  that  the  results  are  not  sensitive  to  this  parameter  in  the 
integration. 

The  integration  was  also  performed  for  two  values  of  (5  and  89°)  correspond¬ 
ing  approximately  to  a  point  charge  and  full  hemispherical  emission. 

When  the  current  at  0  =  90°  due  to  emission  from  the  area  defined  by  a  = 

m 

5°  is  normalized  to  the  emission  area  (Eq.  A-21),  the  numerical  ly  integrated  result  is 


l90.  =0.67^  ,  (A-25) 

which  is  close  enough  to  Eq.  A-20  to  indicate  that  the  numerical  integration  is  reasonable. 

When  the  numerical  ly  integrated  value  for  a  =  89°  is  normalized  to  its  emission 

m 

area,  the  result  is 

l90o  =  1.1 8  ^  .  (A-26) 

Thus,  distributing  the  emitted  charge  q  over  a  hemisphere  increases  the  current  at  6  =  90° 

compared  to  the  current  due  to  a  point  charge  by  a  factor  of  about  1.75. 

9  • 

For  use  with  the  P  formulation,  the  quantity  qv  should  be  identified  with  PA, 

2 

where  A  is  the  emission  area.  For  hemispherical  emission  (A  =  2ttR  ),  the  current  den¬ 
sity  at  9  = 90°  is  then 
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(A  -27) 


H 


90 


,j>(90o)  2ttR 


=  1.18  P 


The  symbol  H  has  been  used  to  denote  the  azimuthal  magnetic  field  which  is  numerically 
equal  to  the  replacement  current  density  at  the  conducting  surface. 

A .4  CAI3LE  COUPLING 


X-ray-induced  cable  si  gna  Is  can  be  i  magined  to  result  from  four  sepa  rate  drivers: 

1.  Knock-on  current  (emitted  electrons  captured  by  the  cable). 

2.  HA  voltage  source. 

3.  CV  current  source. 

4.  Photon-driven  charge  (differential  charge  transfer  between  the  dielec¬ 
tric  and  the  conductors  of  the  cable). 

These  mechanisms  are  defined  and  discussed  in  the  following  sections.  They  are  fairly 
well  understood,  so  the  objective  of  this  study  is  to  derive  quick,  simple  means  for  esti¬ 
mating  the  relative  and,  if  possible,  the  absolute  magnitudes  of  cable  responses  due  to  dif¬ 
ferent  SG  E  M  P  pulses. 

In  this  discussion,  the  word  "cable"  can  mean  the  conductor  of  an  unshielded 
cable,  the  conducting  sheath  of  a  shielded  cable,  or  the  outer  conducting  shield  around  a 
cable  bundle.  The  currents  that  are  estimated  are  the  total  currents  on  the  cables  or 
cable  shields.  It  is  beyond  the  scope  of  the  present  effort  to  atten  pt  to  predict  the  pene- 
t  ration  of  current  through  theshieldof  a  single  cable  or  a  cable  bundle,  or  to  estimate  the 
division  of  the  current  between  separate  wi  res  of  a  cable  bundle. 

In  most  situations,  the  radiation  hardness  analyst  will  not  have  a  precise  def  ini- 
tion  of  the  locations  of  the  cables  in  his  system,  especially  in  the  design  stage  before  the 
final  configuration  is  established.  In  that  case,  he  is  usually  interested  in  worst-case 
values  for  scoping  the  magnitude  of  the  problem  and  for  pinpointing  the  potentially  crit¬ 
ical  items  that  will  need  more  refined  estimates.  Hence,  the  coupling  formulas  used 
herein  assume  worst-case  locations  and  orientations  of  the  cables.  For  example,  for  HA 
coupling,  the  cable  is  assumed  to  be  oriented  parallel  to  the  conducting  surface  and  per¬ 
pendicular  to  the  direction  of  the  local  magnetic  field  Similarly,  for  CV  coupling,  the 
cable  is  assumed  to  be  close  to  the  emitting  surface  where  the  electric  fields  are  largest. 

A. 4.1  Knock-On  Current 

Knock-on  current  is  the  net  flux  of  electrons  captured  or  emitted  by  the  outside 
of  a  cable.  This  current  must  flow  back  to  ground,  mainly  through  the  cable  shield,  but 
some  fraction  of  it  will  couple  into  the  center  conductor  of  the  cable,  depending  on  how 


50 


effective  the  cable  shielding  is.  The  emission  current  from  the  cable  is  estimated  by  the 
simplified  relations  programmed  into  the  71-59  calculator,  ignoring  space-charge-limiting 
from  the  cable  since  cable  sizes  and  distances  to  conducting  structures  are  usually  small 
compared  to  the  width  of  the  space-charge  layer.  Thus,  in  the  present  discussion,  we  will 
be  concerned  only  with  the  amount  of  current  captured  by  a  cable  due  to  emission  from 
adjacent  surfaces. 

For  very  low  fluxes,  when  space-charge-limiting  is  negligible,  the  capture  current 
is  just  the  emission  current  density  from  the  surrounding  surfaces  times  the  appropriate 
projected  area  of  the  cable.  However,  when  there  is  a  significant  space-charge  potential 
barrier  and  the  cable  is  not  immediately  adjacent  to  the  emitting  surface,  the  lower- 
energy  emitted  electrons  cannot  reach  the  cable  because  of  the  potential  barrier.  On  the 
other  hand,  the  higher-energy  electrons  will  be  able  to  reach  the  cable,  and  those  elec¬ 
trons  on  a  direct  line  with  the  cable  will  be  immediately  captured  by  it.  Some  of  the 
nearby  electrons  that  missed  the  cable  on  the  outward  flight  will  be  turned  back  by  the 
potential  barrier  and  can  then  strike  the  cable  on  their  way  back  to  the  emitting  surface. 
Thus,  a  realistic  estimate  of  the  total  capture  current  is  crucially  dependent  on  a  know¬ 
ledge  of  the  potential-barrier  profile  in  the  region  of  the  emitting  surface  (the  space- 
charge  region).  This  problem  has  been  studied  extens  ive  ly  in  the  SGEMP  com  muni  ty ,  bot  h 
for  conducting  cavities  (short-circuited  diodes)  and  for  external  surfaces  (open-circuit 
diodes);  see,  for  example,  References  A-10  -  A-12.  However,  because  of  the  complexity 
of  the  problem,  most  of  the  analytical  studies  have  been  restricted  to  steady-state  condi¬ 
tions.  Obviously,  early  in  a  pulse,  before  the  potential  barrier  has  been  ful  ly  established,  a 
larger  percentage  of  the  emitted  current  can  reach  a  given  cable  position.  Whether  this 
initial  surge  of  current  is  greater  than  the  peak  current  in  a  fully  developed  barrier 
depends  on  the  rate  of  rise  of  the  emission  pulse  and  the  degree  of  space-charge  limiting. 
Reference  A-10  describes  a  complex  computer  code  which  treats  the  time-dependent 
problem,  but  results  are  not  easily  reduced  to  analytic  solutions.  Considerable  simpli¬ 
fication  of  the  problem  has  been  assumed  for  the  calculator  codes  at  present. 

Knock-on  current  is  the  difference  between  (1)  the  current  of  electrons  that  is 
emitted  from  the  illuminated  surface  of  the  satellite  and  is  then  captured  by  the  cable  and 
(2)  the  current  emitted  from  the  outer  surface  of  the  cable  by  the  incident  photons.  For 
worst-case  estimates,  it  is  conservative  to  consider  the  captured  and  emitted  currents 
separately#  since  they  always  produce  cable  currents  of  opposite  signs. 
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For  the  emitted  current,  the  current  per  unit  length  of  cable  can  be  taken  as  the 
circumference  of  the  cable  times  the  larger  of  the  reverse  or  forward  yield  from  the 
cable,  divided  by  the  photon  pulse  width  (FWHM).  These  yields  can  be  obtained  by  the 
routines  described  in  Section  A. 2.  In  all  these  calculations  of  knock-on  current,  the  rele¬ 
vant  cable  radius  corresponds  to  its  outermost  surface,  whether  the  latter  is  a  conductor 
or  a  dielectric. 

For  the  captured  current,  the  magnitude  depends  on  the  height  of  the  potential 
barrie  r  re  I  at  ive  to  the  energy  spect  rum  of  emitted  elect  r  on  s  and  the  location  of  the  cable 
relative  to  the  position  of  the  peak  in  the  potent ial.  However,  for  worst-case  estimates,  it 
can  be  assumed  that  the  cable  is  fairly  close  to  the  emitting  surface  and  that  the  potential 
barrier  is  large  enough  that  practically  all  the  emitted  electrons  are  reflected  back  to  the 
emitting  surface.  Thus,  the  cable  will  capture  electrons  on  both  its  front  and  back  sides. 
Except  for  a  slight  time  delay  for  the  emitted  electrons  to  be  reflected  from  the  potential 
barrier,  the  instantaneous  knock-on  current  due  to  the  captured  electrons  is  just 


Kt)  =4a*)em(t). 


(A -2  8) 


where  4a  A  is  twice  the  projected  area  of  the  cable  and  I  (t)  is  the  instantaneous 
emission  current  density. 

A. 4.2  B  Coupling 

This  coupling  mechanism  is  the  voltage  induced  in  a  conducting  loop  by  the  rate 
of  change  of  the  magnetic  flux  lines  enclosed  by  that  loop. 

As  mentioned  previously,  the  B  fields  on  the  outside  of  a  photo-illuminated  body 
are  proportional  to  P,  the  rate  of  change  of  the  space-charge  dipole  moment.  Thus,  a 
relative  comparison  of  B  coupling  can  be  obtained  by  comparing  the  peak  P's  for  differ¬ 
ent  excitations.  However,  for  system  coupling  studies,  one  needs  absolute  values  of  the 
magnetic  fields,  and  very  often  one  needs  them  inside  conducting  cavities. 

Consider  a  cable  with  the  centerline  above  a  conducting  surface,  h,  and  the 
length  A  defining  the  coupling  area  (A  =hH).  If  the  cable  is  close  to  a  surface,  the 
driving  B  field  is  the  one  due  to  the  returning  skin  current,  Ks  (amp/m),  which  is  caused  by 
the  rate  of  change  of  the  space-charge  dipole  moment  P.  In  that  case,  the  driving  volt¬ 
age  in  the  circuit  is 


dK 

"A  *  M0A  dT 


-  V  KP  p 


(A-29) 
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where  mq  is  the  magnetic  permeability  of  free  space,  K*  is  the  constant  proportionality 
•  ** 
between  P  and  the  skin  current  «s  at  the  point  of  interest  on  the  body,  and  P  would  be 

obtained  as  described  in  Section  A. 3.4. 

« 

For  BA  coupling  inside  a  conducting  cavity,  the  return  currents  on  the  walls  of 
the  cavity  will  be  a  function  of  the  rate  of  change  of  the  space-charge  dipole  moment  P 
in  the  cavity  and  the  portion  of  the  emitted  current  that  surmounts  the  space-charge 
potential  barrier  and  is  transmitted  across  the  cavity.  For  large  cavities,  high  fluences, 
and/or  very  short  pulses,  the  contribution  of  P  to  the  peak  surface  currents  should  exceed 
the  contribution  from  the  transmitted  current.  On  the  other  hand,  for  smal  I  cavities,  low 
fluences,  and/or  long  pulses,  the  surface  currents  will  be  determined  primarily  by  the 
transmitted  current. 

For  BA  coupling  due  to  the  t  ransmitted  current  in  a  conduct ing  cavity,  it  is  con¬ 
servative  for  present  purposes  to  assume  that  the  depth  of  the  cavity  is  small  so  that 
space-charge-limiting  will  be  minimal.  It  is  assumed  that  the  user  will  input  a  character¬ 
istic  dimension  (2R)  of  the  cavity  normal  to  the  incident  beam.  Assuming  that  this  cavity 
is  a  cylindrical  can  of  radius  K  which  is  illuminated  from  one  end,  the  maximum  magnetic 
field,  H,  will  be  adjacent  to  the  cylindrical  wall  and  will  have  a  magnitude 


H  =  'em!  ' 


(A-30) 


where  Jem  is  the  internal  emission  current  density  (amp/m^)  from  the  end  of  the  cylinder. 
The  assumption  of  a  cylindrical  cavity  gives  an  upper-limit  est  imate  of  the  return  current 
densities,  except  for  local  effects,  because  a  cylinder  has  the  smallest  possible  ratio  of 
ci rcumference  to  cross-sectional  (emission)  area.  The  resulting  BA  voltage  is 


d  dJ 

a  a  aii  .  K  em 

BA  =  WQ  AH  =  I'q  A  J  ~dt_ 


( A-3 1) 


The  amount  of  current  that  will  flow  on  the  cable  as  a  result  of  the  BA  voltage 
will  depend  on  its  inductance,  L,  and  its  termination  resistance,  Z  R.  For  a  simple  circuit 
with  L  and  ZR  in  series  (Figure  A-7)  and  driven  by  a  step  voltage  BA,  the  resulting 
current  in  the  cable  is 
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Figure  A-7.  Simplified  equivalent  ci  rcuit  for  BA 
driver  excitation  of  a  cable 


For  a  large  value  of  ZR,  the  exponential  term  damps  out  before  the  driving  voltage  (BA) 
decays  away.  In  that  case,  the  peak  I  is  just  BA/ZR.  However,  for  essentially  short- 
circuit  conditions  (Z  R  -  0),  the  current  I  wil  I  increase  linearly  with  t  ime  for  the  duration 
of  the  voltage  driver  (tp),  so  the  peak  current  is 
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for  a  constant  B. 

It  is  assumed  that  the  user  will  input  his  termination  impedance,  Z^,  or  use  a 
default  value  of  50  ohms.  In  addition,  if  the  short-circuit  condition  applies,  he  must  input 
the  radius  of  his  cable  or  cable  bundle,  a,  in  addition  to  the  previously  inputted  height  of 
the  cable  centerline  above  the  conducting  plane,  h.  The  code  will  then  calculate  the 
inductance  L  from  the  following  formula  for  a  cylindrical  conductor  over  a  ground  plane. 


L 


2* 


£n 


(A-34) 


A .4.3  CV  Coupling 


When  the  voltage  V  between  two  conductors  changes  with  time,  the  resulting 

• 

driver  for  an  equivalent  electrical  circuit  can  be  represented  by  a  current  source  CV, 
where  C  is  the  capacitance  between  the  conductors.  As  mentioned  previously,  for  worst- 
case  estimates  it  is  assumed  that  the  cable  is  close  to  the  emitting  surface  where  the 
electric  fields  are  largest,  either  on  the  exterior  of  a  satellite  or  in  an  interior  cavity.  For 
the  present  purposes,  it  is  conservative  to  ignore  the  decrease  in  the  electric  field  with 
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distance  from  the  emitting  surface,  due  either  to  the  space-charge  cloud  or  to  the  image 
charges  in  the  opposite  wal  I  of  a  cavity. 

Using  the  peak  value  of  the  rate  of  change  of  the  surface  electric  field,  E,  the 
peak  V  is  just  Eh,  where  h  is  the  height  of  the  cable  above  the  emitting  surface.  Using 
the  radius  of  the  cable  a  and  the  cable  length  l,  the  capacitance  C  for  a  cylinder  over  a 
ground  plane  can  be  written  as 


C  = 


2  it  Ic, 


Jtn 


|h  +  /h^  -  a_2  j 


( A — 3  5) 


where  is  the  permittivity  of  free  space.  Note  that  Eqs.  A-34  and  A-35  yield  a  propaga¬ 
tion  velocity  along  the  cable  equal  to  the  speed  of  light,  as  they  should. 

For  the  simple  circuit  shown  in  Figure  A-8,  the  current  I  in  the  cable  due  to  a 
step  current  driver  hEC  is 

l(t)  =  A-j  exp(-s.jt)  +  Aj  expl^t)  +  hEC  ,  (A-36) 


where 


s1  2  =  0.5  ~  t  /(Zr/L)*  -  4c/L 


and  A^  and  A2  are  constants  determined  by  the  boundary  conditions. 


( A  — 3  7) 
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Figure  A-8.  Simplfied  equivalent  circuit  for  CV 
excitation  of  a  cable 
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If  the  duration  of  E  is  long  compared  to  the  damping  times  1/s^  and  l/Sj,  the 
peak  current  in  the  cable  will  be  the  steady-state  value  hEC.  However,  if  R  is  very  large, 
Eq.  A-36  is  approximately 

l(t)  *  hECjl  -  exp[-t(/ZR  C)j|  ,  (A-38) 

so  l(t)  initially  increases  linearly  as 

l(t)  =  hEt/Z  R  ,  (A-39) 

and  the  peak  current  is 

'p  =  hEmax/2R*  (A -40) 

To  scope  the  cable  current  when  the  termination  impedance  is  uncertain,  the 
peak  current  should  be  calculated  both  by  Eq.  A-40  and  from  the  steady-state  value  hEC, 
and  the  larger  value  should  then  be  used. 

A. 4.4  Photon-Oriven  Current 

When  photons  are  incident  on  a  cable  or  a  cable  bundle,  electrons  are  emitted 
from  the  cable  conductors  and  are  trapped  in  the  cable  dielectrics.  This  displacement  of 
charge  causes  a  potential  difference  between  the  cable  core  and  shield,  which  drives  a 
current  between  the  two  conductors.  This  current  is  called  the  photon-driven  current. 

Extensive  work  is  reported  in  the  literature  on  measurements  and  predictions  of 
these  currents  (see,  for  example.  Refs.  A-1  3  -  A-1  5).  The  magnitude  of  the  current,  or  the 
total  charge  transfer  during  the  pulse,  is  a  function  of  the  cable  geometry  and  materials 
and  of  the  energy  spectrum  and  intensity  of  the  emitted  electrons.  An  estimate  of  the 
relative  magnitudes  of  this  current  in  a  given  cable  in  vacuum  due  to  different  photon 
pulses  can  be  made  by  comparing  the  dipole  moment,  D,  of  the  charge  emitted  from  the 
cable  conductors  and  captured  in  the  dielectrics.  This  dipole  moment  is  just  the  product 
of  the  emission  yield  times  the  average  range  of  the  emitted  electrons  in  the  dielectric, 
including  gaps  between  the  cable  conductors  and  the  dielectrics.  The  absolute  magnitude 
of  this  cable  response  can  be  est  imated  by  integrating  the  effect  of  D  around  the  circum¬ 
ference  of  the  cable,  near  both  the  center  and  outer  conductors.  This  integration  has  been 
done  fairly  successfully  for  coaxial  cables  graphically  (Ref.  A-7)  and  with  a  computer  code 
(Ref.  A-15).  However,  both  these  approaches  require  considerable  reduction  to  fit  into 
present  programmable  calculators  such  as  the  TI-59. 
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The  present  expression  employed  on  the  calculator  is 
I  =  2$P£/AT  , 

where  $  is  the  fluence  (cal/m-4),  P  is  the  dipole  moment  per  unit  area  of  the  electrons 
emitted  by  the  shield  into  the  dielectric  material  (C-m/cal),  i  is  the  cable  length  (m), 
and  AT  is  the  x-ray  pulse  width  (sec).  P  is  based  on  the  simple  planar  geometry  value  in 
the  data  base  obtained  with  the  QUICKE2  electron  code.  P  can  be  replaced  by  the  elec¬ 
tron  yield  times  a  gap  width  if  one  exists,  and  if  that  product  is  greater  than  P.  The 
approximation  provides  better  than  order-of-magnitude  estimates,  and  should  scale  prop¬ 
erly  with  x-ray  and  material  parameters  when  compared  with  more  exact  calculations. 

A. 4. 5  Summary  of  Cable  Drivers 

The  cable  coupling  terms  are  summarized  in  Figure  A-9.  The  upper  right  corner 
defines  spectrum  and  geometry  parameters  where  $  is  the  x-ray  flux  (cal/mzAec),  f  indi- 
cates  the  energy  spectrum,  t  is  the  attenuation  material  mass  (g/cm  ),  Z  is  the  atomic 
weight  of  the  electron  emission  material,  lg  is  the  emitted  electron  current  ( amp/m 2),  and 
the  other  input  variables  are  self-explanatory.  L  and  C  are  the  cable  inductance  and 
capacitance  per  unit  length,  respectively.  E  and  B  are  the  electric  and  magnetic  field 
rates. 
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APPENDIX  B 


USER'S  INSTRUCTIONS  FOR 

PROGRAMMABLE  CALCULATOR  CODES  FOR  SCEMP  EFFECTS 

B.1  FOREWORD 

This  appendix  provides  step-by-step  instructions  to  the  SGEMP  analyst  for 
operating  the  calculator  codes  described  in  the  main  body  of  this  report.  The  programs 
are  designed  for  use  on  the  Texas  Instruments  TI-59  calculator  and  do  not  require  a  printer 
or  any  other  special  options.  Codes  and  data  are  read  in  from  magnetic  cards,  and  the 
user  keys  in  selected  parameters  describing  the  x-ray  envi  ronment  and  satellite  structure. 

The  steps  requi  red  by  the  user  are  to  (1)  specify  the  x-ray  spect  rum,  (2)  attenuate 
it  to  the  satellite  location  of  interest,  (3)  compute  SG  E  M  P  excitation  parameters,  (4)  com¬ 
pute  SCEMP  response  parameters,  and  finally  (5)  compute  cable  current  drivers.  Each  of 
these  analysis  functions  is  addressed  by  a  separate  sma II  program  in  the  calculator  which 
is  read  in  at  the  beginning  of  each  step.  Abbreviated  job  descriptions  and  their  related 
calculator  codes  are  listed  in  Table  B-1.  The  codes  are  listed  in  the  logical  flow  of  analy¬ 
sis  from  specification  of  the  x-ray  environment  to  calculation  of  cable  current  drivers. 
The  codes  generally  flow  together  in  that,  where  possible,  data  from  one  program  used  by 
another  does  not  have  to  be  manually  handled  between  steps.  The  second  program  simply 
picks  it  up  from  where  it  is  located  and  proceeds. 

The  user  is  assumed  to  have  elementary  knowledge  of  the  calculator's  operation. 
Some  general  calculator  and  code  characteristics  are  discussed  in  Section  B.2.  The  indi¬ 
vidual  programs  are  described  in  detail  in  Sections  B.3  through  B.7.  A  complete  sample 
problem  is  described  in  Section  B.8. 


Table  B-1.  Calculator  Codes  for  SGEM  P  Analysis  Steps 
(Parentheses  around  a  program  name  indicate  that  it  is  optional) 


Code  Name 

Function 

None  Required 
(TBLSET) 

(BLKBDY) 

ATTEN 

EXCITE 

RESPONSE 

CABLE 

Specify  x-ray  spectrum 

Specify  arbitrary  x-ray  spectrum 

Specify  blackbody  x-ray  spect  rum 

Attenuate  x  rays 

Calculate  electron  yield,  velocity,  and  dipole  moment  in  dielectrics 
Calculate  E&M  response  parameters 

Calculate  cable  current  drivers 

PHBCBDUK)  M0B  fIU0 

I 


B.2 


USFR  INSTRUCTIONS  COMMON  TO  ALL  PROGRAMS 


This  section  describes  procedures  common  to  al  I  of  the  programs  in  this  appendix. 

A  brief  overview  of  the  codes  is  helpful  in  understanding  their  operation.  The. 
calculator  memory  is  divided  into  four  banks  which  can  be  programmed  with  up  to  240 
steps  each.  Memory  can  be  partitioned  to  interchange  code  steps  and  data  storage,  with 
up  to  30  numbers  stored  in  each  bank.  One  insertion  of  a  magnetic  card sets  all  locations 
in  a  bank  of  memory.  As  an  example,  in  the  program  EXCITE,  the  memory  has  been  pro¬ 
grammed  according  to  the  divisions  shown  in  Table  B-2.  The  memory  allocation  shown 
limits  the  program  to  only  240  steps,  but  it  has  the  considerable  advantage  of  permitting 
spectra  to  be  stored  on  one  set  of  cards  and  excitation  parameter  data  on  another  set. 
The  user  may  read  in  a  spectrum  to  bank  2  and  computer  parameters  for  different  satellite 
configurations  by  simply  reading  in  different  data  to  bank  3  and  repeating  the  calculation. 


Table  B-2.  TI-59  Memory  Allocation  for  Computation  of 
Spectrum  Evaluation  Parameters,  Program  EXCITE 


Bank  1 

Program 

Up  to  240  steps 

Bank  2 

Spect  rum 

Up  to  14  point  pai  rs 

Bank  3 

Excitation 
parameter  curves 

Up  to  15  point  pairs 

Bank  4 

Worki  ng  storage 

Magnetic  cards  are  labeled  according  to  information  on  the  top  or  bottom  half  of 
the  magnetic  strip  on  the  back.  The  arrows  in  the  corners  indicate  the  direction  in  which 
the  card  is  to  be  fed  into  the  right-hand  slot  on  the  calculator.  The  card  is  rotated  to  read 
in  the  right-hand  arrows,  not  turned  over  (the  user  can  always  see  the  labeling).  The  upper 
label  on  the  card  pertains  to  the  left-hand  arrow,  the  lower  to  the  right. 

The  calculator  memory  is  partitioned  differently  in  different  sections  of  the  sys¬ 
tem.  Program  and  data  cards  are  marked  in  the  lower  right-hand  corner  in  terms  of  the 
number  of  program  steps  and  data  storage  locations  existing  in  the  memory  at  the  time 
the  card  was  written.  The  default  setting  of  the  calculator  is  479.59,  indicating  480  pro¬ 
gram  steps  and  60  data  storage  locations,  corresponding  to  a  50-50  division  of  the  memory 
(data  storage  location  requires  eight  program-steps'  worth  of  memory).  Cards  written  at 
default  partitioning  are  not  marked.  The  repartitioning  of  the  memory  is  accomplished  by 
pressing  n  Op  17,  where  n  is  the  number  of  decades  of  data  storage  desired.  For  example, 
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many  of  the  programs  and  data  use  the  partitioning  of  239.89,  or  240  program  steps  and  90 
data  locations.  This  is  accomplished  in  memory  by  pressing  9  Op  1  7.  To  examine  the  par¬ 
titioning  at  any  time,  press  Op  16. 

The  calculator  labels  are  the  keys  in  the  top  row  (A-E,  A'-E").  These  are  special- 
function  keys  which  can  be  used  to  perform  tasks  automatical  ly.  For  example,  data  can  be 
stored  in  the  correct  location  automatically  by  the  program  when  the  user  enters  his  num¬ 
ber  and  presses  special-function  key  A,  instead  of  having  to  specify  the  exact  location  in 
memory  (using  a  STO  30  operation,  for  example). 

B.3  X-RAY  SPECTRUM  INPUT 

The  x-ray  spectrum  is  stored  in  the  calculator  in  locations  60-73  (energy  in  keV) 
and  locations  75-88  (differential  fluence  in  units  of  cal/cm^AeV).  These  values  can  be 
entered  manually,  or  with  optional  routines  described  below,  or  they  may  be  already  in 
existence  on  a  magnetic  card.  Toperform  manual  ent  ry,  values  for  energies  and  fluences 
are  stored  by  hand,  with  the  user  keeping  track  of  the  storage  location.  To  input  a  spec¬ 
trum  without  having  to  simultaneously  keep  track  of  the  location  within  the  calculator, 
refer  to  the  instructions  for  the  program  TBLSET  later  in  this  section.  To  specify  a  black- 
body  spectrum  automatically,  see  the  instructions  for  program  BIKBDY.  To  read  in  a 
spectrum  already  on  a  data  card,  partition  the  calculator  to  the  same  partitioning  it  had 
when  the  card  was  written,  press  CLR,  and  insert  the  card.  The  calculator  will  display 
-2,"  indicating  that  memory  bank  2  has  been  loaded. 

Spectra  must  have  at  least  three  energy  bins,  and  must  have  a  zero  following  the 
last  entry  (the  zero  should  be  in  location  74  for  spectra  having  the  maximum  of  14  energy 
bins).  The  present  data  base  for  photoemission  excitation  parameters  is  limited  to  the 
range  1  to  1000  keV,  so  x  rays  should  be  defined  in  this  range  only.  Normalization  to  unit 
incident  fluence  is  accomplished  automatically  by  the  E  XCITE  program. 

To  save  a  spect  rum  on  a  card,  press  CLR,  WRITE,  and  insert  a  card  to  be  writ¬ 
ten.  The  card  may  be  subsequently  read  to  specify  a  spectrum  by  simply  pressing  CLR  and 
inserting  the  card,  as  described  above.  It  is  suggested  that  the  calculator  partitioning  at 
the  time  of  writing  the  card  be  placed  on  the  card  to  prompt  the  user  attempting  to  read 
in  that  spect  rum. 

The  recipes  for  different  spectrum  input  methods  are  given  in  the  following. 
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B.3.1 


Enter  a  value  for  energy  into  the  calculator  and  press  STO  n,  where  n  is  the 
appropriate  storage  location,  60  <  n  <  73.  Enter  a  value  for  differential  fluence  and  press 
ST O  m,  where  m  is  the  appropriate  location,  75  <  m  <  88.  Bin  1  energy  is  in  location  60, 
differential  fluence  in  75,  bin  2  energy  in  61,  etc.  Observe  general  rules  for  spectra  men¬ 
tioned  in  the  beginning  of  Section  B.3. 

B.3.2  Optional  Spectrum  Input  Program:  TBLSET 

A  routine  is  available  to  expedite  table  entries  and  checkouts.  The  routine  is 
designed  for  generalized  table  entry  at  arbitrary  locations  in  memory.  In  the  present  case, 
a  string  of  energies  must  be  placed  into  memory  beginning  at  location  60,  and  the  corre¬ 
sponding  differential  fluences  beginning  at  location  75. 

To  perform  this  operation,  press  9  Op  17  to  partition  memory,  then  CL  R,  then 
read  in  the  TBLSET  code  card.  Press  RST,  then  R/S.  The  calculator  now  displays  1.,  indi¬ 
cating  that  the  calculator  is  ready  to  receive  the  energy  value  for  energy  bin  1.  Energies 
may  now  be  input  by  keying  them  in  and  pressing  R/S.  The  index  of  the  next  bin  then 
appears.  To  indicate  the  last  bin,  enter  a  zero  and  the  calculator  will  automatically  skip 
to  location  75  and  display  1.,  indicating  that  the  calculator  is  ready  to  receive  differential 
fluence  for  bin  1.  At  this  point,  begin  entering  differential  fluences  and  pressing  R/S. 
Stop  after  the  last  bin. 

The  values  may  be  displayed  in  succession  by  pressing^  and  then  R/S,  R/S,  etc. 
Displays  beginning  at  any  location  can  be  commenced  by  entering  the  desired  location  and 
pressing  A' .  The  calculator  will  automatically  begin  displaying  at  location  60  if  no  entry 
is  made  into  A'  . 

A  summary  of  TBLSET  operations  is  given  in  Table  B-3. 

Table  B -3.  Summary  of  Input/Output  Routing  TBLSET  Inst  ructions 


Labels:  A' 

Starting  loca¬ 
tion  for  display 
(default  = 
value  in  B) 

B' 

c'  \r 

Starting  loca¬ 
tion  to  multiply, 

R/S,  multiply 
by  R/S 

E  ' 

B 

C  D 

E 

Starting  loca- 

Display  va lues 

tion  for  storage 

beginning 

(default  =  60) 

at  A' 
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Table  B-3.  Summary  of  Input/Output  Routing  TBLSET  Instructions  (Cont.) 


Press  9  Op  17. 

Press  C  L  R. 

Read  in  program  TBLSET. 

Press  R ST,  R/S. 

Enter  beginning  storage  location  in  B  (default  value  =  60). 

Enter  energy,  press  R/S,  energy,  R/S,  etc. 

Terminate  string  with  0.,  R/S. 

Enter  differential  fluence,  press  R/S,  etc. 

To  display  values,  press _E_,  R/S,  R/S,  .  .  . 

To  display  values  beginning  at  a  location  other  than  B,  enter  the  location  in  A' 
after  RST,  then  press  _E_,  R/S,  R/S,  .  .  . 

To  save  the  spectrum  on  a  card,  press  CL  R ,  2_,  WRITE,  and  insert  a  magnetic 
card. 


B.3.3  Optional  Blackbody  Spectrum  Generator  Program:  BLKBDY 
Setup 

Partition  the  calculator  to  240  program  steps  and  90  data  locations. 

Press  9  Op  17.  239.89  will  be  displayed. 

Press  C  L  R. 

Insert  program  card  (BLKBOYO.  1.  will  be  displayed. 

Input 

A.  Number  of  bins  —  i.e.,  the  number  of  (energy,  differential  fluence)  pairs  desired; 
a  minimum  of  three  bins  and  a  maximum  of  14  bins  can  be  specified. 

Enter  the  number  of  bins  desired  into  the  calculator  and  press  A. 

B.  Blackbody  temperature  (keV). 

Enter  the  blackbody  temperature  into  the  calculator  and  press  jl. 

C.  Normalization  fluence,  $  (optional). 

Units  as  user's  choice. 

The  blackbody  spectrum  is  normalized  so  that  the  total  fluence  equals  4>.  (The 
default  is  4>  =  1.) 

Enter  the  total  fluence  $  into  the  calculator  and  press^_. 

E  xecution 


Press  E  to  execute. 


L  ocat  ion 

Energies  (keV)  60-73 

Differential  fluences  (cal/cm^/keV)  75-88 

To  save  the  spect  rum  on  a  magnet  i  c  card,  press  C  L  K,  2_,  WRITE,  and  insert  a 
card. 

Cautions 

A.  For  blackbody  temperatures  T  greater  than  80  keV,  the  spectrum  obtained  with 
BLKRDY  could  give  erroneous  results  when  used  in  EXCITE.  Thi s  problem  arises 
because  the  blackbody  spectrum  for  T  >  80  will  contain  energies  greater  than 
1000  keV  and  the  excitat  ion  parameterdatabaseusedin  EXCITE  is  valid  on  ly  for 
energies  from  1  to  1000  keV. 

B.  If  BLKBDY  is  executed  twice  in  a  row  and  the  number  of  bins  is  changed,  the 
blackbody  temperature  and  normalization  constant  must  a  Iso  be  re-entered. 

B.4  ATTENUATION  PROGRAM:  ATEN 

Setup 

Partition  calculator  to  240  program  steps  and  90  data  locations. 

Press  9  Op  17.  239.89  wil  I  be  displayed 

Press  C  L  R. 

Insert  program  card  (AT  EN ).  1.  will  be  displayed 

Input 

A.  Energy  spectrum  input  --  see  Section  B.3. 

B.  Attenuation  material  data. 

1.  X-ray  cross-section  data:  Data  necessary  for  calculation  of  the  cross  sec¬ 
tion  is  stored  in  locations  21-29.  Cards  are  available  for  low-Z  (Z  x  29)  and 
high-Z  (Z  >  29)  materials,  where  Z  is  the  atomic  number.  To  read  in  the 
card,  press  CL  R  and  insert  the  card  appropriate  to  the  material  of  interest 
(either  high-  or  low-Z).  The  calculator  will  display  4.,  indicating  that  the 
data  has  been  stored  in  memory  bank  4. 

2 

2.  Mass  per  unit  area,  x  (g/cm  ),  of  the  attenuation  material:  To  read  in  x, 
enter  the  value  for  x  into  the  calculator  and  press  _A_. 

3.  Attenuation  material  atomic  number,  Z:  To  read  in  Z,  enter  its  value  into 
the  calculator  and  press  C. 
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E  xecut ion 

The  calculator  is  ready  to  calculate  the  total  incident  fluence,  the  attenuated  spec¬ 
trum,  and  the  total  attenuated  fluence. 

A.  Total  incident  fluence,  4>|n  (cal/cm^)  (optional).  If  is  desired,  press  E'.  At 
the  termination  of  the  calculation,  4>|n  is  displayed.  This  optional  step  integrates 
the  incident  spect rum  without  modifying  it. 

R.  Attenuated  spectrum.  Press  When  the  display  is  fixed,  the  at  tenuated  spec¬ 
trum  has  been  stored  in  the  appropriate  data  registers.  This  step  alone  simply 
stores  the  values  and  does  not  print  out  the  integrated  fluence. 

C.  Total  attenuated  fluence,  1  n  (cal/cm^)  (optional).  If  <I>^ttnis  desired,  press 
E  "  after  R  above  is  accomplished.  At  the  termination  of  the  calculation,  is 

displayed.  This  step  is  done  after  R  above  if  it  is  done  at  all. 

As  an  example,  if  the  user  has  entered  his  spectrum  and  wants  to  know  the  incident 
and  attenuated  fluences,  he  presses  E E ,  and  then  E  "  again.  If  he  only  wants  the 
attenuated  fluence,  he  presses  ^and  then  E  If  he  intends  to  execute  EXCITE  for  his 
next  step,  the  integrated  fluence  is  computed  automatically  and  step  is  not 
necessary  in  AT  EN. 

Attenuation  Through  Several  Materials 

A  given  incident  spectrum  can  be  attenuated  through  any  number  of  materials. 
For  example,  for  two  materials,  material  1  and  material  2:  First  input  the  incident 
energy  spectrum  and  attenuation  material  data  for  material  1.  Execute  the  program. 
The  attenuated  spectrum  from  material  1  is  now  the  incident  spectrum  of  material  2 
and  is  already  stored  in  the  proper  data  locations.  Therefore,  enter  the  appropriate 
attenuation  material  data  for  material  2.  When  the  program  is  executed,  the  result¬ 
ing  spectrum  has  been  attenuated  through  both  material  1  and  material  2.  This  can  be 
repeated  for  any  number  of  materials. 

R.5  EXCITATION  PARAMETER  PROGRAM:  EXCITE 

Setup 

Partition  calculator  memory  to  240  maximum  program  steps  with  00  data 
locations: 

Press  9  Op  1 7 
Press  C  L  R . 

Insert  program  card  (called  EXCITE). 

Calculator  wil  I  display  1.,  indicating  that  bank  1  has  been  programmed. 
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Excitation  Parameter  Data  Input 

Read  selected  excitation  parameter  data  into  locations  30-59  (bank  3): 

Press  CL R. 

Insert  card. 

Calculator  displays  3.,  indicating  that  bank  3  has  been  programmed. 

E  xecution 

The  calculator  is  now  ready  to  compute  the  selected  parameter  for  the  spectrum 
chosen. 

Press  RST. 

Press  R/S. 

Execution  commences  and  numbers  flash  on  the  display  for  each  energy  bin  in  the 
order  below. 

e.  i**1  bin  energy  (keV) 

f  h  y 

4>0e  j  iin  differential  fluence  (cal/cm  /keV) 

T.  $oe  •  i  cumulative  sum  of  fluences  in  each  bin  (cal/cm2) 

i 

When  the  display  becomes  fixed,  the  value  is  total  fluence: 

Nbins 

]£  ^Oei  Aei  (cal/cm2), 

i  =1 

where  Nbins  is  the  total  number  of  energy  bins.  This  is  the  total  fluence  incident 
at  the  emission  surface  of  the  excitation  parameter  configuration. 

Press  R/S. 

Calculator  displays  excitation  parameter  P  (per  unit  fluence  on  the  emitting 
material).  P  represents  any  excitation  parameter. 

Press  R/S. 


Warning:  EXCITE  gives  the  photoelectron  yield,  for  example,  per  unit  fluence  at  the 
emission  location.  The  x-ray  shinethrough  fraction  must  then  be  multiplied  in  to  obtain 
the  total  yield  for  x-ray  fluence  on  the  outside  of  the  attenuating  material. 
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Calculator  displays  ratio  Pprevjous'/Pcurrent'  t*ie  rat'°  °f  l^e  previously  calcu¬ 
lated  excitation  parameter  to  the  current  one.  Thus,  if  the  velocity  is  being 
computed,  for  example,  the  yield  (for  velocity  calculation)  is  computed  first, 
then  a  new  card  read  in  for  yield/velocity  and  the  process  repeated.  The  third 
static  display  yields  the  ratio  yield/fyield/velocity)  =  velocity. 

B.6  RESPONSE  PARAMETER  PROGRAM:  RESPONSE 

Setup 

Response  is  stored  on  two  magnetic  cards  and  occupies  banks  1-3  in  the  calcu¬ 
lator.  Banks  1  and  2  are  stored  on  one  card  and  bank  3  on  the  second  card.  The 
calculator  must  be  partitioned  to  560  program  steps  and  50  data  locations. 

Press  5  Op  17.  559.49  wil  I  be  displayed. 

Press  CL  R. 

Insert  program  card  RESPONSE,  bank  1.  1.  will  be  displayed. 

Press  CL  R. 

Insert  program  card  RESPONSE,  bank  2.  2.  will  be  displayed. 

Press  CL  R. 

-  Insert  program  card  RESPONSE,  bank  3.  3.  will  be  displayed. 

Input 

The  input  parameters  must  be  manually  stored  into  the  following  storage 
locations. 

Parameter  Storage  Location 


*R 

Pulse  rise  t  ime  (sec) 

10 

At 

FWHM,  pulse  full  width  at  half  maximum  (sec) 

12 

4> 

Total  fluence  (cal/m  ) 

13 

Y 

Y ield  (C/cal) 

16 

v 

ave 

Average  normal  velocity  (mAec) 

17 

For  example,  if  tR  =  10  nsec, 

Press 

1  E  E  V-  8. 

Press  STO  10. 

Proceed  similarly  for  At ,  <•>,  Y,  and  vgve. 
E  xecution 

To  execute  RESPONSE: 

Press  _PST« 

Press  R/S. 
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Output 

The  calculator  display  will  become  fixed  when  all  parameters  have  been  calcu¬ 
lated.  RESPONSE  calculates  values  for  the  output  parameters  and  stores  the- 
values  in  the  following  locations. 

Parameter  Storage  Location 


^  em 

2 

Peak  emission  current  (amp/m  ) 

15 

• 

P 

Peak  time  rate  of  change  of  the 
space-charge  dipole  moment  (amp/m) 

20 

P 

• 

Peak  t  i me  rate  of  change  of  P  (amp/m-sec) 

21 

• 

E 

Peak  time  rate  of  change  of  electric  field 
(V/m-sec) 

22 

E 

Peak  electric  field  (V/m) 

23 

lm1 1 

lm2  i 

18 

■  Time  of  peak  P  (whichever  is  non-zero) 

19 

lm3 

Time  of  peak  E,  time  for  P  =  0  (sec) 

26 

*lf  the  value  for  E  is  zero,  the  problem  was  computed  to  be  non-space-charge-limited,  for 
which  no  formula  for  E  has  provided  in  the  programs.  The  value  for  E  is  simply 
I  E  At/e,  however,  in  that  case. 

For  example,  the  value  for  P  can  be  displayed  in  the  calculator  by  pressing  RC  L 
2U 

RESPONSE  can  be  executed  for  other  input  parameters  by  simply  changing  the 
ones  being  varied  and  pressing  RST,  R/S.  It  is  not  necessary  to  reenter  param¬ 
eters  held  constant. 

B.7  CARLE  CURRENT  PROGRAM:  CARLE 
Setup 

The  program  CABLE  is  stored  on  two  sides  of  one  magnetic  card  and  occupies 
banks  1  and  2  in  the  calculator.  The  calculator  must  be  partitioned  to  560  pro¬ 
gram  steps  and  50  data  locations. 

Press  5  Op  17.  559.49  will  be  displayed. 

Press  CLR. 

Insert  program  card  CABL  E,  bank  1  1.  will  be  displayed. 

Press  CLR. 

Insert  program  card  CABL  E,  bank  2  2.  will  be  displayed. 
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Input 

A.  Response  quantities.  The  necessary  response  quantities  can  be  entered  by  run¬ 
ning  the  program  RESPONSE  or  they  can  be  manually  stored  into  the  calculator 
memory. 

1.  Inputs  from  RESPONSE:  CABLE  has  been  designed  to  use  as  input  results 
obtained  by  executing  the  program  RESPONSE.  All  necessary  response 
parameters  are  calculated  by  RESPONSE  and  stored  in  the  correct  locations 
for  use  in  C  ABL  E. 

2.  The  necessary  response  quantities  can  also  be  manually  entered  into  the  cal¬ 
culator.  These  response  parameters  must  be  stored  in  the  following  loca¬ 
tions.  None  of  these  parameters  needs  to  be  specified  if  RESPONSE  has 
been  executed  previously. 

Parameter  Location 


eo 

Permittivity  of  free  space  =  8.854  x  10"12 

00 

1 R 

Pulse  rise  time  (sec) 

10 

At 

FWHM,  pulse  full  width  at  half  maximum  (sec) 

12 

*  em 

2 

Peak  emission  current  (amp/m  ) 

15 

• 

P 

Peak  time  rate  of  change  of  the 

space-charge  dipole  moment  (amp/m) 

20 

P 

• 

Peak  time  rate  of  change  of  P  (amp/m-sec) 

21 

• 

E 

Peak  time  rate  of  change  of  electric  field 
(V/m-sec) 

22 

E 

Peak  electric  field  (V/m) 

23 

Inputs 

The  cable  inputs  are  the  cable  standoff,  cable  length,  table  radius,  termination 
impedance,  and  cavity  diameter.  Not  all  of  these  parameters  are  necessary 
inputs  in  calculating  the  various  currents.  Consult  Section  B.7  to  determine 
which  inputs  are  required  for  the  current(s)  of  interest. 

The  cable  inputs  are  entered  as  follows. 

Cable  standoff,  h  (m) 

Enter  value  for  h. 

Press  A'  . 
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a 


Cable  length,  i  (m) 

Enter  value  for  A. 

Press  C' . 

Termination  impedance,  ZR  (ohm) 

Enter  value  for  ZR  (value  must  be  zero). 
Press  O'- 

Cable  radius,  a  (m) 

Enter  value  for  a. 

Press  E". 

Cavity  diameter,  2R  (m) 

Enter  value  for  2R . 

Press  E. 


E  xecution 

Listed  below  are  execution  instructions  for  the  various  currents  and  inputs 
needed  for  each  current  (see  paragraph  B  under  "Input"). 

A.  Surface  skin  current  density  (numerically  equal  to  tangential  magnetic  f  ieldH  )  , 

<P 

H,  (amp/m) 

<P 

Press  A.  The  value  forH,  is  displayed  and  also  is  stored  in  location  26. 

.  ~  * 

B.  BA  currents,  I  •  (amp) 

•  ext 

E  xterior  BA  current,  I* 
n 

Cable  inputs  needed  —  h,  A ,  Z  R ,  a 

Press  B.  The  value  for  I.,  is  displayed  and  stored  in  location  27. 

—  t  •* 

Interior  BA  current,  l.*n 
H 

Cable  inputs  needed  --  h,  l,  2R,  Z  R,  a 

Press  B'.  The  value  for  I  ^nt  is  displayed  and  stored  in  location  28. 

C.  CV  current,  I  •  (amp) 

Cable  inputs  needed  —  h,  l,  ZR,  a 

Press  JC.  The  value  for  I  •  is  displayed  and  stored  in  location  29. 

D.  Knock-on  current,  l|<noc|<  (amp) 

Cable  inputs  needed  —  a 

PressJD.  The  value  for  I^q^  is  displayed  and  stored  in  location  30. 
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E.  Direct  drive  current,  Iqq  (amp) 

Cable  inputs  needed  --  l 

Enter  the  value  for  the  dipole  moment  P  (C-cm/cal),  obtained  from 

EXCITE. 

Press  RST. 

Press  R/S.  The  value  for  Iq^  is  displayed  and  stored  in  location  31. 

B.8  SAMPLE  PROBLEM 

X  rays  characterized  by  an  8-keV  blackbody  spectrum  are  incident  on  a  3-mil- 
thick  aluminum  box.  The  total  incident  fluence,  4>.|n,  is  5  cal/m2  and  the  pulse  is  charac¬ 
terized  by  a  12-nsec  rise  time  and  a  FWHM  of  20  nsec.  Calculate  the  peak  time  rates  of 
change  of  the  electric  field  and  space-charge  dipole  moment.  Calculate  the  cable  cur¬ 
rents  I  •  and  I  •  for  the  configurations  given  in  Figure  B-1. 

H  E 


RE-03279 


Figure  B-1.  Sample  problem  configuration 
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This  problem  can  be  solved  quickly  and  easily  on  the  TI-59  calculator  using  the 
programs  described  in  the  previous  sections.  The  method  of  solution  involves  the  following 
five  steps. 

1.  Cenerate  an  8-keV  blackbody  spectrum  using  the  program  BL  KHl)Y. 

2.  Attenuate  the  spectrum  through  3  mils  of  aluminum  using  the  program 
ATEN. 

3.  Calculate  excitation  parameters  (electron  emission  yield  and  normal  veloc¬ 
ity)  using  the  program  EXCITE. 

4.  Calculate  the  response  parameters  (peak  time  rates  of  change  of  the  electric 
field  and  space-charge  dipole  moment)  using  the  program  RESPONSE. 

5.  Couple  these  response  parameters  into  the  system  and  calculate  cable  cur¬ 
rents  using  the  program  CABLE. 

Step-by-step  instructions  used  in  solving  this  sample  problem  are  given  below  for 
each  of  the  f  ive  programs  mentioned  above.  The  f  ive  programs  are  to  be  executed  sequen¬ 
tially  since  required  input  for  one  program  is  calculated  and  stored  by  the  previous  pro¬ 
gram.  In  general,  if  a  mistake  is  made  during  input  of  parameters,  it  is  safe  tostart  over 
at  the  beginning  of  the  ca Icul at iona I  s  ep  after  magnetic  cards  are  read  in  for  that  step. 

B.8.1  BLKBDY  (see  Section  R.3) 

Generate  an  8-keV  blackbody  spectrum  with  total  fluence$^n=  5  cal/m^.  An 
energy  spectrum  is  specified  by  a  curve  of  differential  fluence,  $  (cal/cm  -keV)  versus 
energy,  E  (keV),  and  thus  consists  of  (Ej,$j)  pairs.  In  general,  the  greater  the  number  of 
pairs  (bins),  the  more  accurately  a  given  spectrum  can  be  represented.  A  maximum  of  14 
(E  j,<I>  j)  pairs  can  be  used  in  these  programs,  and  if  carefully  chosen,  is  enough  to  accurately 
represent  most  spectra  of  interest.  However,  for  this  sample  problem,  only  five  (Ej,$j) 
pai  rs,  wil  I  be  calculated,  to  shorten  the  t  ime  requi  red  to  run  through  this  example. 
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E  nter 

Press 

Insert 

Display 

Comments 

CLR 

0 

9 

Op  17 

239.89 

Partition  calculator 

CLR 

0 

BLKBOY 

1. 

• 

Read  in  magnetic  card 

5 

A 

5. 

5  ( E j , 4> j )  pairs  specified 

8 

B 

1. 

8-keV  spectrum  specified 

5 

C 

5. 

Total  fluence  <t>^.  =  5  cal/mz 

E 

.0130935239 

Execute;  calculates  and  stores  (E  j,<t> j )  pairs 

• 

The  calculator  occasionally  has  difficulty  reading  cards,  especially  if  they  were  written 

by  a  different  calculator.  A 

flashing  display  at  the  end  of  this  operation  means  that  a  read 

problem  occurred. 

and  the  card  should  be  inserted  again  after  pressing  CLR. 

The  generated  spectrum  which  is  stored  in  locations  60-88  can  be  examined  as 

follows. 

E  nter 

Press 

Insert 

Display 

Comments 

RCL  60 

3.2 

> 

0> 

s - ' 

r— 

LU 

RCL  61 

24. 

E2  (keV) 

RCL  62 

00 

• 

E3  (keV) 

RCL  63 

72. 

E4  (keV) 

RCL  64 

96. 

E5 (keV) 

RCL  65 

0. 

0.  stored  in  location  following  final  energy 

RCL  75 

.0130935239 

4'1  (cal/m2-keV) 

RCL  76 

.1423463343 

4>2  (cal/m2-keV) 

RCL  77 

.0540071936 

4^  (cal/m2-keV) 

RCL  78 

.0090535248 

4.  (cal/m2-keV) 

4 

RCL  79 

.0010683155 

4,.  (cal/m2-keV) 

B.8.2  AT EN  (see  Section  8.4) 

Attenuate  the  8-keV  blackbody  spectrum  through  3  mils  of  aluminum.  The  den- 

1 

sity  of  aluminum  is  p  =  2.7  g/cm  ,  so  i,  the  mass  per  unit  area,  is  calculated  to  be 
t  ss  (3  x  1 0'3  inch)(2.54  cm/fnch)(2.7  g/cm^), 
x  =  0.021  g/cm2. 
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Enter 


Press  Insert 


Display  Comments 


CL  R 

0 

ATEN 

1. 

Read  in  magnetic  card 

CL  R 

0 

LOW  Z 

4. 

Read  in  mag.  card;  Ail  (Z=13),  low-Z  material 

.021 

A 

.021 

t  =  0.021  g/cm2 

13 

C 

13 

Atomic  number  for  Ai,  l  =13 

E  ' 

5. 

Total  incident  fluence,  <t>|n  =  5  cal/m2 

E 

0 

Execute;  the  attenuated  spectrum  is 
calculated  and  stored 

E' 

4.621826108 

Total  attenuated  fluence,  tn  =  4.62182610H 

cal/tn2  T 

The  attenuated  spectrum  is  stored  in  locations  60-88  and  can  be  examined  as 

follows. 

Enter 

Press 

Insert 

Display 

Comments 

RCL  60 

3.2 

E1  (keV) 

RCL  61 

24. 

E2  (keV) 

RCL  62 

48. 

E3 (keV) 

RCL  63 

72. 

rr, 

a'—'. 

ID 

< 

RCL  64 

96. 

E5  (keV) 

RCL  65 

0. 

0.  stored  in  location  following  final  energy 

RCL  75 

0.00000001 8 

**ttn  <cal/m2-keV) 

RCL  76 

.1378648063 

4>^ttn  (cal/m2-keV) 

RCL  77 

.0537916668 

*fttn  (cal/m2-keV) 

RCL  78 

.0090428046 

4>|ttn  (cal/m2-keV) 

RCL  79 

.0010677817 

$*ttn  (cal/m2-keV) 

B.8.3  EXCITE  (see  Section  B.5) 

Calculate  the  forward  electron  emission  yield  and  average  normal  velocity  in  the 
following  manner. 
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1 


Enter  Press 

Insert 

Display 

Comments 

CLR 

0 

EXCITE 

1. 

Head  in  magnetic  card 

CL  R 

0 

YIELD  (for 
velocity 
calculation) 

3. 

Read  in  magnetic  card;  yield  data  base 
(for  velocity  calculation) 

RST 

3. 

R/S 

4.621826108 

Total  fluence  of  spectrum  in  locations  60-88 

R/S 

« 

0000000745 

Yield  per  unit  total  fluence 

R/S 

Store  yield  (ignore  display) 

CLR 

0 

YIELD 

Read  in  mag.  card,  yield/velocity,  data  base 

VELOCITY 

RST 

3. 

K/S 

4.621826108 

Total  fluence  of  spectrum  in  locations  60-88 

R/S 

1 

.5368727-15 

Yield/velocity  per  unit  total  fluence 

R/S 

48447820.7 

Yield/(yield/velocity)  =  velocity 

CLR 

0 

FORWARD 

YIELD 

ALUMINUM 

3. 

Read  in  magnetic  card,  forward  yield  data 
base  for  aluminum 

RST 

3. 

R/S 

4.621826108 

Total  fluence  of  spectrum  in  locations  60-88 

R/S 

0000000745 

Yield  per  unit  total  fluence 

Results: 

Average  normal  velocity  =  48447820.7  m/sec  =  4.84  x  1()7  m/sec 


Yield  per  unit  incident  fluence  = 


_ yield _ 

attenuated  fluence 


attenuated  fluence 
incident  fluence 


yield _ 

incident  fluence 


=  (.0000000745) 


4.621826108 

5.0 


=  6.89  x  10'8  C/cal 
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F 


B.8.4  RESPONSE  (see  Section  B.6) 

Calculate  the  peak  time  rate  of  change  of  the  electric  field  E  and  the  peak  time 
rate  of  change  of  the  space-charge  dipole  moment  P. 


E  nter 

Press  Insert 

Display 

Comments 

CLK 

0 

5 

Op  17 

559.49 

CLK 

0 

RESPONSE 
(Bank  1) 

1. 

Read  in  magnetic  card. 

CLK 

0 

RESPONSE 
(Bank  2) 

2. 

Read  in  magnetic  card. 

CLK 

0 

RESPONSE 
(Bank  3) 

3. 

Read  in  magnetic  card. 

12  EE 

12  00 

9 

*/- 

12-09 

STO  10 

1.2-08 

Pulse  rise  time,  t^  (sec). 

20 

EE 

20  00 

9 

V- 

20-09 

STO  12 

2.0-08 

FWHM,  At  (sec). 

5 

STO  13 

5.00 

2 

Q  Total  incident  fluence  (cal/m  ) 

6.89 

EE 

6.89  00 

8 

*/- 

6.89-08 

STO  16 

6.89-08 

2 

Yield  per  incident  fluence,  Y  (C/cal-m  ). 

4.84 

EE 

4.84  00 

7 

STO  17 

4.84  07 

Average  normal  velocity,  v  (m/sec). 

K  ST 

4.84  07 

R/S 

8.3369  08 

E  xecute. 

The  calculated  response  parameters  can  be  examined  as  follows. 
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Enter 


Press  Insert 


Display  Comments 


CLK 

0 

RCL  15 

17.225 

2 

Peak  emission  current,  lem  (amp/rn  ). 

RCL  20 

9.307368187 

Peak  time  rate  of  change  of  space-charge 
dipole  moment  (am p/m). 

RCL  21 

833690000. 

• 

Peak  t  ime  rate  of  change  of  P,  P  (amp/m-sec). 

RCL  22 

1.9454484  12 

• 

Peak  time  rate  of  change  of  electric  field,  E 

(V/m-sec). 

RCL  23 

0. 

For  tm3  >  At ,  program  sets  E  =0. 

RCL  18 

0. 

• 

Since  value  is  zero,  use  tm2  as  t  ime  of  peak  P  . 

RCL  19 

.0000000228 

• 

Time  of  peak  P,  tm2  =  22.8  nsec. 

RCL  26 

.0000000324 

• 

Timeofpeak  E  =  32.4  nsec. 

B.8.5 

CABLE  (see  Sect  ion  B. 7) 

Calculate  the  cable  currents  1  •  and  1  •  for  the  conf  igu rat  ion  in  Figure  B-1. 

H  t 

Enter 

Press 

Insert 

Display 

Comments 

CLR 

0 

CABLE 
(Rank  1) 

1. 

Read  in  magnetic  card. 

CLR 

0 

CABLE 
(Bank  2) 

2. 

Read  in  magnetic  card. 

.02 

A' 

0.02 

Cable  standoff,  h  (m). 

.2 

C' 

2.  -01 

Cable  lengt  h,  i  (m). 

0 

nr 

0.  00 

Termination  impedance,  Z^  (ohm). 

.002 

E" 

2.  -03 

Cable  radius,  a  (m). 

.2 

E 

2.  -01 

Cavity  diameter,  2R  (m). 

B' 

3.6157638-02 

l£  (amp). 

.05 

A' 

5.  -02 

Cable  standoff,  h  (m). 

.15 

C' 

1.5  -01 

Cable  length,  l  (m). 

.001 

E' 

1.  -03 

Cable  radius,  a  (m). 

C 

1.762642-01 

1  •  (amp) 

Results: 


I ,•  =  0.04  amp 
I  •  =  0.18  amp 

B.8.6  Bin  Selection 

As  mentioned  earlier,  a  large  number  of  bins  [i.e.,  (Ej,$j)  pairs]  may  be  required 
to  accurately  represent  a  given  spectrum.  In  this  sample  problem,  only  five  bins  were 
specified  to  save  time.  For  comparison,  the  same  sample  problem  was  run  specifying  14 
bins  in  BLKBDY.  Results  of  these  two  calculations  are  given  below. 


Result 

5  bins 

14  bins 

Average  velocity  (m/sec) 

4.84  x  107 

3.72  x  107 

Yield  (C/cal) 

6.89  x  10-8 

1.12  x  10"7 

'  em 

17.225 

28. 

• 

P 

9.31 

7.92 

P 

8.34  x  108 

7.62  x  108 

• 

E 

1.95  x  1012 

3.16  x  1012 

■h 

0.04 

0.06 

'e 

0.18 

0.29 

There  are  definite  differences  between  the  results  of  these 

emphasized  that  the  user  must  select  bins  with  care. 

two  calculations,  and  it 

APPENDIX  C 

PHOTOEMISSION  EXCITATION  PARAMETER  DATA  BASE 


This  appendix  contains  graphs  of  the  photoemisson  excitation  parameter  data 
base  employed  in  the  calculator  programs.  The  materials  and  excitation  parameters  for 
which  data  are  provided  on  magnetic  cards  are  listed  in  Table  C-1.  The  average  normal 
electron  velocity  is  only  a  weak  function  of  the  emission  material  and  direction,  and  in 
fact,  the  forward-directed  emission  data  for  aluminum  is  representative,  within  a  factor 
of  2,  of  every  material  for  both  forward  and  backward  directions.  The  dipole  moment  data 
are  provided  for  representative  cable  shield  materials.  Extension  of  the  data  base  to  other 
materials  and  parameters  can  be  performed  straightforwardly  as  the  need  arises.  Rele¬ 
vant  quantities  such  as  dose  and  electron  range  can  be  included,  for  example. 

Results  for  electron  yield,  average  normal  velocity,  and  dipole  moment  in  dielec¬ 
tric  materials  are  plotted  in  Figures  C-1  through  C-7  as  a  function  of  in-and-out  x-ray 
energy  over  the  range  1  to  1000  keV.  The  curves  were  obtained  from  the  QUICKE2  elec¬ 
tron  emission  code  for  monoenergetic  x-ray  spectra,  and  are  fit  to  within  a  factor  of  2 
over  the  energy  range  using  power-law  interpolation  in  the  calculator.  As  many  as  14 
point  pai  rs  in  both  the  x-ray  spectrum  and  the  data  base  are  interpolated  and  convolved  by 
the  EXCITE  program  (see  Appendix  B)  to  obtain  the  total  quantity  for  a  given  spectrum. 
The  units  of  electron  yield  and  dipole  moment  are  presented  as  coulombs  per  calorie  and 
coulomb-meters  per  calorie.  The  yield  is  actually  coulombs  per  unit  area  per  calorie  per 
unit  area,  which  reduced  to  coulombs  per  calorie.  The  dipole  moment  can  be  thought  of 
similarly. 
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Table  C-1.  Materials  and  Excitation  Parameters  Available 
on  Magnetic  Cards  for  the  T  1-59  Calculator  Program  EXCITE 


Photo-E lect  ron 

Y  ield 

Average 

Dipole  Moment  in 

Low-Z  Material 

Material 

Forward 

Reverse 

Velocity 

Forward  Reverse 

Carbon 

X 

X 

|| 
i/i  r  • 

«  §  <N 

Aluminum 

X 

X 

^  o 

x  E  o 

X  X 

Copper 

X 

X 

Qj  +* 

ta  u 

>  0) 

Z.  </> 

X  X 

Silver 

X 

X 

,£  <D  * 

E  l£ 

X  X 

Gold 

X 

X 

§  ■-  -= 
i-c.tr 
-  c  J 

Mylar 

X 

X 

E  *  _ 

3  -O  ** 

~  »-  n3 

5': 

Fiberglass 

X 

X 

o  !  * 

^  C  w 

m  ^  r 

Quartz 

X 

X 

—  'JZ 

U 


C_J 

o 


RE-03280 


Figure  C-1.  Reverse-emit  ted  photoelectron  yield  versus 
monoenergetic  x-ray  energy  for  represent  at  ive  elements 
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i  WMM» 


■»  — 1  ' 


YIELD  (C/cal) 


PHOTON  ENERGY  (keV) 


Figure  C-2.  Reverse-emitted  photo¬ 
electron  yield  versus  monoenerget ic 
x-ray  energy  for  representative  sat¬ 
ellite  materials 
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Figure  C-3.  Electron  average  nor¬ 
mal  velocity  versus  monoenergetic 
x-ray  energy  used  to  represent  all 
elements  and  satellite  materials  for 
forward  and  reverse  electron 
emission 
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1000 


Figure  C-4.  Forward-emitted  photo 
electron  yield  versus  monoenergetic 
x-ray  energy  for  elements 


Figure  C-5.  Forward-emitted  photo 
electron  yield  versus  monoenergetic 
x-ray  energy  for  representative 
satellite  materials 
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4  cy  ATTN:  TITL 

Defense  Technical  Information  Center 
12  cy  ATTN:  DD 

Field  Command 
Defense  Nuclear  Agency 

ATTN;  FCP,  J.  T.  McDaniel 

ATTN:  FCLMC,  H.  R.  Putnam 

ATTN:  FCT 0 

ATTN:  FCTT 

ATTN:  FCTT,  W.  Surma 

ATTN:  FCTP 

ATTN:  FCTXE 

Field  Command 
Defense  Nuclear  Agency 
Livermore  Branch 
ATTN:  FCPRL 

Field  Command 
Defense  Nuclear  Agency 
Los  Alamos  Branch 
ATTN:  FCPRA 

Field  Command  Test  Directorate 
Test  Construction  Division 
Defense  Nuclear  Agency 
ATTN:  FCTC 

Interservice  Nuclear  Weapons  School 
ATTN :  TTV 

Joint  Chiefs  of  Staff 

ATTN:  0-5,  Nuclear  Division 
ATTN:  C3S,  Evaluation  Office 

Joint  Strat  Tgt  Planning  Staff 
ATTN:  JLTW-2 
ATTN:  JLA 

National  Conmunications  System 
ATTN:  NCS-TS 

Under  Secretary  of  Defense  for  Rsch  A  Engrg 
ATTN:  Strategic  &  Space  Sys  (OS) 
ATTN:  AE 


DEPARTMENT  OF  THE  ARMY 

BMD  Advanced  Technology  Center 
Department  of  the  Army 
ATTN:  ATC-0 

BMD  Systems  Command 
Department  of  the  Army 
ATTN:  BDMSC-H 

Deputy  Chief  of  Staff  for  Rsch  Dev  8.  Acq 
Department  of  the  Army 
ATTN:  DAMA-CSS-N 

Electronics  Tech  &  Devices  Lab 
U.S.  Army  Electronics  R&D  Command 
ATTN:  DRSEL 

Harry  Diamond  Laboratories 
Department  of  the  Army 

ATTN:  DELHD-N-RBC ,  R.  Gilbert 
ATTN:  DELHD-I-JL 

U.S.  Army  Communications  Sys  Agency 
ATTN:  CCM-AD-LB 

U.S.  Army  Foreign  Science  &  Tech  Ctr 
ATTN:  DRXST-IS-1 

U.S.  Army  Missile  Command 
ATTN:  RSIC 

DEPARTMENT  Of  THE  NAVY 

Naval  Research  Laboratory 

ATTN:  Code  7550,  J.  Davis 

ATTN:  Code  6 701 

ATTN:  Code  6707,  K.  Whitney 

Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
ATTN:  NSP 

DEPARTMENT  OF  THE  AJR  FORCE 

Air  Force  Geophysics  Laboratory 
ATTN:  PH,  C.  Pike 

Air  Force  Weapons  Laboratory 
Air  Force  Systems  Command 
ATTN:  SUL 
ATTN:  NT 

ATTN:  NTSW,  Co)  Bloemker 
2  cy  ATTN:  NTYC 

Air  University  Library 
Department  of  the  Air  Force 
ATTN:  AIIL/LSE 
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DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 

Ballistic  Missile  Office 
Air  Force  Systems  Command 
ATTN:  ENMG 
ATTN:  SYDT 
ATTN:  ENSN 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
Department  of  the  Air  Force 
ATTN:  AFRDQ1 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  YKD 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN :  YNV 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN:  ESR,  E.  Burke 

Strategic  Air  Command 
Department  of  the  Air  Force 
ATTN:  XPFS 

ATTN:  NRI-STINFO,  Library 

OTHER  GOVERNMENT  AGENCIES 

Central  Intel  1 igence  Agency 

ATTN:  OSWR/STD/MTB ,  A.  Padgett 

Department  of  Commerce 
National  Oceanic  &  Atmospheric  Admin 
ATTN:  F.  Fehsenfeld 

NASA 

Lewis  Research  Center 
ATTN:  Library 
ATTN:  N.  Stevens 
ATTN:  C.  Purvis 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Lawrence  Livermore  National  Lab 

ATTN:  Technical  Info  Dept,  Library 

Los  Alamos  National  Laboratory 
ATTN :  MS  364 

Sandia  National  Laboratories 
Livermore  Laboratory 
ATTN:  T.  Dellin 

Sandia  National  Lab 
ATTN:  3141 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  V.  Josephson 
ATTN:  Library 

Avco  Research  &  Systems  Group 
ATTN:  Library  A830 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
Beers  Associates,  Inc 


ATTN 

B.  Beers 

Computer  Sciences  Corp 

ATTN 

A.  Schiff 

Dikewood  Corporation 

ATTN 

Technical  Library 

Dikewood  Corporation 

ATTN 

K.  Lee 

EG&G  Wash  Analytical  Svcs  Ctr, 

ATTN 

Library 

Ford  Aerospace  &  Communications 

ATTN- 

A.  Lewis 

General  Electric  Co 

ATTN 

J.  Peden 

Hughes  Aircraft  Co 

ATTN 

Technical  Library 

Hughes  Aircraft  Co 

ATTN 

E.  Smith 

ATTN 

W.  Scott 

ATTN 

A.  Narevsky 

Institute  for  Defense  Analyses 

ATTN 

Classified  Library 

IRT  Corp 

ATTN 

N.  Rudie 

ATTN 

B.  Williams 

ATTN 

Library 

JAYCOR 

ATTN 

Library 

4  cy  ATTN 

E.  Wenaas 

4  cy  ATTN 

A.  Woods 

4  cy  ATTN 

R.  Leadon 

4  cy  ATTN 

B.  Nichols 

4  cy  ATTN 

K.  Brannon 

JAYCOR 

ATTN 

R.  Sullivan 

JAYCOR 

ATTN 

R.  Poll 

Johns  Hopkins  University 

ATTN 

P.  Partridge 

Kaman  Sciences  Corp 

ATTN 

Library 

ATTN 

W,  Rich 

ATTN 

N.  Beauchamp 

ATTN 

D.  Osborn 

Kaman  Tempo 

ATTN 

W.  McNamara 

ATTN 

DAS  I AC 

Lockheed  Missiles  A  Space  Co,  1 

ATTN 

Dept  85-85 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

McDonnell  Douglas  Corp 
ATTN:  R.  Kloster 

McDonnell  Douglas  Corp 

ATTN:  S.  Schneider 

Mission  Research  Corp 

ATTN:  C.  Longmire 
ATTN:  M.  Scheibe 
ATTN:  R.  Stettner 

Mission  Research  Corp 
ATTN:  B.  Goplen 

Mission  Research  Corp,  San  Diego 
ATTN:  Library 
ATTN:  V.  Van  Lint 

Mission  Research  Corporation 
ATTN:  U.  Ware 

Pacific-Sierra  Research  Corp 
ATTN:  L.  Schlessinger 
ATTN:  H.  Brode 

R  &  D  Associates 

ATTN:  S.  Siegel 

ATTN:  Technical  Information  Center 
ATTN:  P.  Haas 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Rockwell  International  Corp 
ATTN:  Library 

Science  Applications,  Inc 
ATTN:  W.  Chadsey 

Science  Applications,  Inc 
ATTN:  D.  Woodward 
ATTN:  1C.  Sites 

Spire  Corp 

ATTN:  R.  Little 

SRI  International 

ATTN:  Library 

Systems,  Science  &  Software,  Inc 
ATTN:  A.  Wilson 
ATTN:  Library 

Teledyne  Brown  Engineering 
ATTN:  I.  Ford 

TRW  Defense  &  Space  Sys  Group 
ATTN:  D.  Clement 

ATTN:  Technical  Information  Center 
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